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Abstract. A new method is described for constructing a generalized solution for 
stochastic differential equations. The method is based on the Cameron-Martin ver- 
sion of the Wiener Chaos expansion and provides a unified framework for the study 
■ of ordinary and partial differential equations driven by finite- or infinite-dimensional 

noise with either adapted or anticipating input. Existence, uniqueness, regularity, 

t"**- ■ and probabilistic representation of this Wiener Chaos solution is established for 

a large class of equations. A number of examples are presented to illustrate the 

i_i general constructions. A detailed analysis is presented for the various forms of the 

passive scalar equation and for the first-order ltd stochastic partial differential equa- 

p I ■ tion. Applications to nonlinear filtering if diffusion processes and to the stochastic 

—J ' Navier-Stokes equation are also discussed. 
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1. Introduction 
Consider a stochastic evolution equation 

(1.1) du(t) = (Au(t) + f(t))dt + (Mu{t) + g{t))dW{t), 

where A and M are differential operators, and W is a noise process on a stochastic basis F = 
(Q., J 7 , {^(}i>o,P). Traditionally, this equation is studied under the following assumptions: 

(i) The operator A is elliptic, the order of the operator M is at most half the order of 
.A, and a special parabolicity condition holds. 

(ii) The functions / and g are predictable with respect to the filtration {J~t\t>0i an d 
the initial condition is .T-o-measurable. 

(hi) The noise process W is sufficiently regular. 

Under these assumptions, there exists a unique predictable solution u of (jl.lj) so that 
u 6 £2(0 x (0) T);H) for T > and a suitable function space H (see, for example, Chapter 
3 of 42 ). Moreover, there are examples showing that the parabolicity condition and the 
regularity of noise are necessary to have a square integrable solution of 

The objective of the current paper is to study stochastic differential equations of the type 
(jl.lj) without making the above assumptions (i)— (iii). We show that, with a suitable def- 
inition of the solution, solvability of the stochastic equation is essentially equivalent to 
solvability of a deterministic evolution equation dv = (Av + <p)dt for certain functions ip; 
the operator A does not even have to be elliptic. 

Generalized solutions have been introduced and studied for stochastic differential equations, 
both ordinary and with partial derivatives, and definitions of such solutions relied on various 
forms of the Wiener Chaos decomposition. For stochastic ordinary differential equations, 
Krylov and Veretennikov |2()j used multiple Wiener integral expansion to study Ito diffusions 
with non-smooth coefficients, and more recently, LeJan and Raimond [221 use d a similar 
approach in the construction of stochastic flows. Various versions of the Wiener chaos 
appear in a number of papers on nonlinear filtering and related topics [21 ESI ESI E23 HH1 
etc.] The book by Holden et al. presents a systematic approach to the stochastic 

differential equations based on the white noise theory. See also ^01) [ID] and the references 
therein. 

For stochastic partial differential equations, most existing constructions of the generalized 
solution rely on various modifications of the Fourier transform in the infinite-dimensional 
Wiener Chaos space L2(W) = L2{p,, Jjf , P). The two main modifications are known as the 
S-transform |10j and the Hermite transform ^2j. The key elements in the development of 
the theory are the spaces of the test functions and the corresponding distributions. Several 
constructions of these spaces were suggested by Hida ^0]) Kondratiev JZJj an d Nualart 
and Rozovskii [HHl- Both S- and Hermite transforms establish a bijection between the 
space of generalized random elements and a suitable space of analytic functions. Using 
the S-transform, Mikulevicius and Rozovskii [33] studied stochastic parabolic equations 
with non-smooth coefficients, while Nualart and Rozovskii [3H] and Potthoff et. al 40_ 
constructed generalized solutions for the equations driven by space-time white noise in 
more than one spacial dimension. Many other types of equations have been studied, and 
the book ^2] provides a good overview of literature the corresponding results. 

In this paper, generalized solutions of ([1.1[) are defined in the spaces that are even larger 
than Hida or Kondratiev distribution. The Wiener Chaos space is a separable Hilbert 
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space with a Cameron-Martin basis P|. The elements of the space with a finite Fourier 
series expansion provide the natural collection of test functions 2?(L2(W)), an analog of the 
space Z>(R ) of smooth compactly supported functions on R d . The corresponding space 
of distributions 2?'(L2(W)) is the collection of generalized random elements represented by 
formal Fourier series. A generalized solution u = u(t, x) of is constructed as an element 
of D'(L2(W)) so that the generalized Fourier coefficients satisfy a system of deterministic 
evolution equations, known as the propagator. If the equation is linear the propagator is 
a lower-triangular system. We call this solution a Wiener Chaos solution. 

The propagator was first introduced by Mikulevicius and Rozovskii in (32]) and further stud- 
ied in [2Hj, as a numerical tool for solving the nonlinear filtering problem. The propagator 
can also be derived for certain nonlinear equations; in particular, it was used in |311 1341 I35j 
to study the stochastic Navier-Stokes equation. 

The propagator approach to defining the solution of (jLljl has two advantages over the S- 
transform approach. First, the resulting construction is more general: there are equations 
for which the Wiener Chaos solution is not in the domain of the S-transform. Indeed, it is 
shown in Section^] that, for certain initial conditions, equation du = u x dWt has a Wiener 
Chaos solution for which the S-transform is not defined. On the other hand, by Theorem 
18.11 below, if the generalized solution of (jl.ip can be defined using the S-transform, then 
this solution is also a Wiener Chaos solution. Second, there is no problem of inversion: the 
propagator provides a direct approach to studying the properties of Wiener Chaos solution 
and computing both the sample trajectories and statistical moments. 

Let us emphasize also the following important features of the Wiener Chaos approach: 

• The Wiener Chaos solution is a strong solution in the probabilistic sense, that is, 
it is uniquely determined by the coefficients, free terms, initial condition, and the 
Wiener process. 

• The solution exists under minimal regularity conditions on the coefficients in the 
stochastic part of the equation and no special measurability restriction on the input. 

• The Wiener Chaos solution often serves as a convenient first step in the investigation 
of the traditional solutions or solutions in weighted stochastic Sobolev spaces that 
are much smaller then the spaces of Hida or Kondratiev distributions. 

To better understand the connection between the Wiener Chaos solution and other notions 
of the solution, recall that, traditionally, by a solution of a stochastic equation we understand 
a random process or field satisfying the equation for almost all elementary outcomes. This 
solution can be either strong or weak in the probabilistic sense. 

Probabilistically strong solution is constructed on a prescribed probability space with a 
specific noise process. Existence of strong solutions requires certain regularity of the coeffi- 
cients and the noise in the equation. The tools for constructing strong solutions often come 
from the theory of the corresponding deterministic equations. 

Probabilistically weak solution includes not only the solution process but also the stochastic 
basis and the noise process. This freedom to choose the probability space and the noise 
process makes the conditions for existence of weak solutions less restrictive than the similar 
conditions for strong solutions. Weak solutions can be obtained either by considering the 
corresponding martingale problem or by constructing a suitable Hunt process using the 
theory of the Dirichlet forms. 



4 



S. V. LOTOTSKY AND B. L. ROZOVSKII 



There exist equations that have neither weak nor strong solutions in the traditional sense. 
An example is the bi-linear stochastic heat equation driven by a multiplicative space-time 
white noise in two or more spatial dimensions: the irregular nature of the noise prevents 
the existence of a random field that would satisfy the equation for individual elementary 
outcomes. For such equations, the solution must be defined as a generalized random element 
satisfying the equation after the randomness has been averaged out. 

White noise theory provides one approach for constructing these generalized solutions. The 
approach is similar to the Fourier integral method for deterministic equations. The white 
noise solution is constructed on a special white noise probability space by inverting an 
integral transform; the special structure of the probability space is essential to carry out 
the inversion. We can therefore say that the white noise solution extends the notion of 
the probabilistically weak solution. Still, this extension is not a true generalization: when 
the equation satisfies the necessary regularity conditions, the connection between the white 
noise and the traditional weak solution is often not clear. 

The Wiener chaos approach provides the means for constructing a generalized solution 
on a prescribed probability space. The Wiener Chaos solution is a formal Fourier series 
in the corresponding Cameron-Martin basis. The coefficients in the series are uniquely 
determined by the equation via the propagator system. This representation provides a 
convenient way for computing numerically the solution and its statistical moments. As a 
result, the Wiener Chaos solution extends the notion of the probabilistically strong solution. 
Unlike the white noise approach, this is a bona fide extension: when the equation satisfies 
the necessary regularity conditions, the Wiener Chaos solution coincides with the traditional 
strong solution. 

After the general discussion of the Wiener Chaos space in Sections |1] and [51 the Wiener 
Chaos solution for equation Ql.l|) and the main properties of the solution are studied in 
Section Several examples illustrate how the Wiener Chaos solution provides a uniform 
treatment of various types of equations: traditional parabolic, non-parabolic, and antici- 
pating. In particular, for equations with non-predictable input, the Wiener Chaos solution 
corresponds to the Skorohod integral interpretation of the equation. The initial solution 
space 2?'(W) is too large to provide much of interesting information about the solution. 
Accordingly, Section discusses various weighted Wiener Chaos spaces. These weighted 
spaces provide the necessary connection between the Wiener Chaos, white noise, and tra- 
ditional solutions. This connection is studied in Section El In Section EH the Wiener Chaos 
solution is constructed for degenerate linear parabolic equations and new regularity results 
are obtained for the solution. Probabilistic representation of the Wiener Chaos solution 
is studied in Section 1101 where a Feynmann-Kac type formula is derived. Sections 111! 1121 
1131 and El discuss the applications of the general results to particular equations: the Zakai 
filtering equation, the stochastic transport equation, the stochastic Navier-Stokes equation, 
and a first-order Ito SPDE. 

The following notation will be in force throughout the paper: A is the Laplace operator, 
Di = d/dxi, i = 1, . . . ,d, and summation over the repeated indices is assumed. The space 
of continuous functions is denoted by C, and H^, 7 6 K, is the Sobolev space 




where / is the Fourier transform of /. 
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Below is a summary of the Hilbert space theory of linear stochastic parabolic equations. 
The details can be found in the books [IJ and [12]; see also For a Hilbert space X, 
(•, -)x and || • || x denote the inner product and the norm in X. 

Definition 2.1. The triple (V,H,V') of Hilbert spaces is called normal if and only if 

(1) V <— > H V and both embeddings V ► H and H V' are dense and continuous; 

(2) The space V is the dual of V relative to the inner product in H ; 

(3) There exists a constant C > so that \(h, v)h\ < C||v||y ||/&||y for all v E V and 
h£H. 

For example, the Sobolev spaces (H e 2 +1 (R d ), H f 2 (R d ), H^ J (R d )), 7 > 0, £ E R, form a 
normal triple. 

Denote by (v 1 , v), v' E V, v E V, the duality between V and V relative to the inner product 
in H. The properties of the normal triple imply that |(u',i>)| <C||u||y||u / ||v 7 ) and, if v' E H 
and v £ V, then (v',v) = (v',v)h', 

Let F = (fi, T, {J"t}t>0i IP) be a stochastic basis with the usual assumptions. In particular, 
the sigma-algebras J- and J-q are P-complete, and the filtration {J~t}t>o 1S right-continuous; 
for details, see |231 Definition 1. 1.1]. We assume that F is rich enough to carry a collection 
Wk = Wk(t), k > 1, t > of independent standard Wiener processes. 

Given a normal triple (V,H,V') and a family of linear bounded operators A(t) : V — * V' , 
Aik(t) :V^H,t£ [0,T], consider the following equation: 



(2.1) u(t) = n + / (^u(s) + /(s))ds + / (M k u(s) + g k (s))dw k (s), < t < T, 

Jo Jo 

where T < 00 is fixed and non-random and the summation convention is in force. 
Assume that, for all v E V, 

(2.2) ^2\\M k (t)v\\ 2 H <oo, tE [0,T]. 

k>l 

The input data uq, f, and g k are chosen so that 

(2.3) E f \\u f H + [ T \\f(t)\\v>dt + V [ T \\9k(t)f H dt I < 00, 

\ Jo k>i Jo J 

uq is .Fo-nieasurable, and the processes f,g k are ^-adapted, that is, f(t) and each g k {t) are 
J-t-measurable for each t > 0. 

Definition 2.2. ^4n Tt-adapted process u E L2(F; L2((0, T); V)) is called a traditional, or 
square-integrable, solution of equation h2.1)) if, for every v E V , there exists a measurable 
sub-set £1' of fl with P(O') = 1, so that, the equality 

(2.4) (u(t),v) H = (u ,v) H + f {Au(s) + f(s),v)ds + J2(M k u(s)+g k (s),v) H dw k (i 

Jo 1 , 

holds on Q' for all < t < T . 



k>l 



Existence and uniqueness of the traditional solution for (J2.1|) can be established when the 
equation is parabolic. 
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Definition 2.3. Equation \2.1\) is called strongly parabolic if there exists a positive num- 
ber e and a real number Co so that, for all v G V and t G [0, T], 

(2.5) 2(A(t)v,v) +J2\\M(t) k 

v \\h + £ \\ v \\v ^ ColMltf- 

k>l 

Equation \2. 1\) is called weakly parabolic (or degenerate parabolic) if condition \2. 5\) holds 
with e = 0. 

Theorem 2.4. If \2. 3\) and \2. 5\) hold, then there exists a unique traditional solution of 
\2.1\) . The solution process u is an element of the space 

L 2 (F; L 2 ((0, T); V)) f| L 2 (F; C((0, T),H)) 

and satisfies 

E ( sup \\u{t)f H + [ \\u(t)\\ldt 

\0<t<T JO 



(2.6) 



<C(C Q ,S,T)M\\\ Uo \\ 2 H + f T \\f(t)f vl dt + Y, l' T \\9k 
\ Jo k>i Jo 



2 H dt 



Proof. This follows, for example, from Theorem 3.1.4 in |42j . □ 

A somewhat different solvability result holds for weakly parabolic equations |42l Section 
3.2]. 

As an application of Theorem 12.41 consider equation 

du(t, x) = (aij(t, x)DiDju(t, x) + bi(t, x)Diu(t, x) + c(t, x)u(t, x) + f(t, x))dt 
+ (a ik (t,x)Diu(t,x) + v k (t,x)u(t,x) + g k (t,x))dw k (t) 
with < t < T, x G M. d , and initial condition u(0,x) = uo(x). Assume that 

(CL1) The functions bounded and Lipschitz continuous, the functions bi, c, Oi k , 

and v are bounded measurable. 
(CL2) There exists a positive number e > so that 

{2a ij (x) - <T ik (x)a jk (x))yiyj > £\y\ 2 , x,y G R d , t G [0,T]. 

(CL3) There exists a positive number K so that, for all x £ R rf , X)fe>i l^fc( x )| 2 ^ K- 
(CL4) The initial condition no £ L 2 (f}; L 2 (R d )) is .T-o-measurable, the processes / G 
L 2 (U x [0,7];ff^ 1 (R d )) and # fc G L 2 (£3 x [0, T]; L 2 QR d )) are ^-adapted, and 

Efc>l/o r]E ll^llL 2 (Md)( t ) dt < °°- 

Theorem 2.5. Under assumptions (CL1)-(CL4), equation \2. 7U /ias a unique traditional 
solution 

u G L 2 (F; L 2 ((0, T); F, 1 (M d ))) f| L 2 (F; C((0, T), L 2 (M d ))), 
and i/ie solution satisfies 



{21 



io ' * k>1 Jo 
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Proof. Apply Theorem 12.41 in the normal triple 

(H^(R d ),L 2 (R d ),H^ 1 (R d )); condition jZHJ) in this case is equivalent to assumption (CL2). 
The details of the proof are in |42[ Section 4.1]. □ 

Condition Q2.5|) essentially means that the deterministic part of the equation dominates the 
stochastic part. Accordingly, there are two main ways to violate (|2.5j) : 

(1) The order of the operator Ai is more than half the order of the operator A. Equation 
du = u x dw(t) is an example. 

(2) The value of *£ k \\M k (t)v\\ 2 H is too large. This value can be either finite, as in 
equation du(t,x) = u xx (t,x)dt + 5u x (t, x)dw(t) or infinite, as in equation 

(2.9) du(t,x) = Au(t,x)dt + a k (x)udw k , a k - CONS in L 2 (R d ), d>2. 

Indeed, it is shown in j3H] that, for equation (|2.9[) . we have 

£l|Mk(tM!r = °° 
fc>i 

in every Sobolev space H 1 . 

Without condition (|2.5|) . analysis of equation (|2.1|) requires new technical tools and a dif- 
ferent notion of solution. The white noise theory provides one possible collection of such 
tools. 

3. White Noise Solutions of Stochastic Parabolic Equations 

The central part of the white noise theory is the mathematical model for the derivative of the 
Brownian motion. In particular, the ltd integral J * f(s)dw(s) is replaced with the integral 
Jo f( s ) oW(s)ds, where W is the white noise process and o is the Wick product. The white 
noise formulation is very different from the Hilbert space approach of the previous section, 
and requires several new constructions. The book 10 is a general reference about the white 
noise theory, while ^2] presents the white noise analysis of stochastic partial differential 
equations. Below is the summary of the main definitions and results. 

Denote by S = S(Mr) the Schwartz space of rapidly decreasing functions and by S' = S'(M. e ), 
the Schwartz space of tempered distributions. For the properties of the spaces S and S' see 

Definition 3.1. The white noise probability space is the triple 

where B(S') is the Borel sigma-algebra of subsets of S' , and [x is the normalized Gaussian 
measure on B{S'). 

The measure \x is characterized by the property 

Js> 

where (cj,</j), u> G 5', (p £ S, is the duality between S and S' . Existence of this measure 
follows from the Bochner-Minlos theorem |121 Appendix A]. 
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Let {rjk, k > 1} be the Hermite basis in L2QR ), consisting of the normalized eigenfunctions 
of the operator 

(3.1) A = -A + |x| 2 , xER £ . 
Each r]k is an element of S |121 Section 2.2]. 

Consider the collection of multi-indices 

Ji = {a = (aj, i>l), aj G {0, 1,2, . . .}, y^Oj < ooj. 

i 

The set J7i is countable, and, for every a £ J, only finitely many of a, are not equal to 
zero. For a G J7i, write a! = ]X a^! and define 

(3.2) = -i=TT J ff a( (( W ,T«)), wG5', 

where (•, •) is the duality between S and S' , and 

(3.3) H B (t) = (-l)V 2 / 2 ^e-^ 

is n th Hermite polynomial. In particular, H\{t) = 1, H\(t) = t, H2(t) = t 2 — 1. If, for 
example, a = (0, 2, 0, 1, 3, 0, 0, . . .) has three non-zero entries, then 

Sa(w) = 21 31 • 

Theorem 3.2. T/ie collection a G i7i} is an orthonormal basis in £2(8). 

Proof. This is a version of the classical result of Cameron and Martin [2] . In this particular 
form, the result is stated and proved in |121 Theorem 2.2.3]. □ 

By Theorem 13. 21 every element ip of £2(8) is represented as a Fourier series <p = ^ a ^ a , 
where <p a = J s , tp(u)£a(u)dn, and |M|| 2(§) = EaeJi l^*' 2 - 

For a G i7i and g £ R, we write 

(2N)" Q = JJ(2j) 9aj - 

i 

Definition 3.3. For p G [0, 1] and q>0, 

(1) i/ie space (S) Pi „ is the collection of elements ip from £2(8) so that 

M% q = E(« ! ) P (2Nn^| 2 <oo; 

(2) the space (S)- p - q is the closure 0/^2(8) relative to the norm 
(3-4) IM| 2 - A - g = ^(a!)-"(2N)-^|^| 2 ; 

(3) the space (S) p is the projective limit of (S) pq as q changes over all non-negative 
integers; 

(4) the space (<S)-p is the inductive limit of (<S)_p_ 9 as g changes over all non-negative 
integers. 



It follows that 
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• For each p £ [0, 1] and q > 0, ((<S) P)9 , L<2(§), (<S)_ P) _ g ) is a normal triple of Hilbert 
spaces. 

• The space (<S) P is a Frechet space with topology generated by the countable family 
of norms || • \\ P) n, n = 0, 1, 2, . . ., and p £ (S) p if and only if <p G (<S) P)9 for every 
q > 0. 

• The space (<S)_ P is the dual of (<S) P and tp £ (<S)- P if and only if (p G (S)- p - q for 
some q > 0. Every element ip from (S) p is identified with a formal sum X^aeJi <p a £a 
so that (j3.4|) holds for some g > 0. 

• For < p< 1, 

c (5) p c (5) c L 2 (S) c (5)_ c (5)_ p c 
with all inclusions strict. 

The spaces (<S)o and (5)i are known as the spaces of Hida and Kondratiev test functions. 
The spaces (5)-o and (<S)_i are known as the spaces of Hida and Kondratiev distributions. 
Sometimes, the spaces (<S) p and (<S)_ P , < p < 1, go under the name of Kondratiev test 
functions and Kondratiev distributions, respectively. 

Let h £ S and hk = J R e h(x)r}k(x)dx. Since the asymptotics of n th eigenvalue of the operator 
A in (|3,1|) is n 1 ^ HU Chapter 21] and A k h £ S for every positive integer k, it follows that 

(3.5) ^2\h k \ 2 k q < oo 

k>i 

for every g G R. 

For a E Si and /i^ as above, write ft" = IIj(^j) aj '> an d define the stochastic exponential 

(3-6) S(h) = £ 

at/a 

Lemma 3.4. TTte stochastic exponential £ = £(/i), h £ S, has the following properties: 

• G (<S) P? < p < 1; 

• For every q > 0, i/iere exists a 5 > so i/iai S(h) G (<S)i i9 as Zon<? as ^fc>i l^fcl 2 < ^- 

Proof. Both properties are verified by direct calculation |12l Chapter 2]. □ 

Definition 3.5. The S-transform Sp(h) of an element ip = ^2 a£ j p a £a from (<S)- P is the 
number 

(3.7) Stp(h) = V —^=<Pa, 

where h = J2k>i ^kVk £ <S and h a = Y\j(hj) aj . 

The definition implies that if (p G (<S)_ P) _ g for some g > 0, then Sp(h) = (p,£(h)), where 
(•,•) is the duality between (<S) P)9 and (5)_ Pi _ g for suitable q. Therefore, if p < 1, then 
S(p(h) is well-defined for all h £ S, and, if p = 1, the S<p(h) is well-defined for /i with 
sufficiently small Z^O^) norm. To give a complete characterization of the S-transform, one 
additional construction is necessary. 

Let U p , < p < 1, be the collection of mappings i 7 from S to the complex numbers so that 
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1. For every h±, hi £ S, the function F(h\+zh2) is an analytic function of the complex 
variable z. 

2. There exist positive numbers K\,Ki and an integer number n so that, for all h £ S 
and all complex number z, 

\F(zh)\ < i^exp <^K 2 \\K n h\\^ d) \z\^ . 

For p = 1, let U l be the collection of mappings F from S to the complex numbers so that 

1'. There exist e > and a positive integer n so that, for all h\,h,2 £ S with 
||A n /ii||£ 2 ( K ^ < e, the function of a complex variable z i— > F(h\ + /12-2) is analytic 
at zero, and 

2'. There exists a positive number K so that, for all h £ S with || ^ n h\\ L2 ^e^ < e, 
|F(/i)| < K. "" 

Two mappings F, G with properties 1' and 2' are identified with the same element of U 1 if 
F = G on an open neighborhood of zero in S. 

The following result holds. 

Theorem 3.6. For every p £ [0, 1], i/ie S-transform is a bijection from (<S)-p to W. 

In other words, for every ip £ (S)-„, the S-transform Sip is an element of U p , and, for every 
F £ U p , there exists a unique 93 S (<S)_ P so that fifyj = F. This result is proved in ^01 when 
p = 0, and in .IT] when p = 1. 

Definition 3.7. For ip and tp from (S)- p , p E [0, 1], the Wick product ipo ip is the unique 
element of (<S)_ P whose S-transform is Sip ■ Sip. 

If S~ l is the inverse S-transform, then 

ip o ip = S~ 1 (Sp ■ Sip), 

Note that, by Theorem lH.6| the Wick product is well defined, because the space U p , p £ [0, 1] 
is closed under the point-wise multiplication. Theorem 13.61 also ensures the correctness of 
the following definition of the white noise. 

Definition 3.8. The white noise W on R e is the unique element of (S)o whose S transform 
satisfies SW(h) = h. 

Remark 3.9. If g £ L p (S), p > 1, then g £ (<S)_o 03 Corollary 2.3.8], and the Fourier 
transform 

g{h) = / exp (y/^l{uj, h)) g(u)dp,(u) 
JS' 

is defined. Direct calculations ^2j Section 2.9] show that, for those g, 

Sg(y/^lh)=g(h)e^ M2 ^). 

As a result, the Wick product can be interpreted as a convolution on the infinite- dimensional 
space (S)- p . 

In the study of stochastic parabolic equations, i = d + 1 so that the generic point from 
M. d+1 is written as (t,x), t £ R, x £ W 1 . As was mentioned earlier, the terms of the type 
fdW(t) become foWdt. The precise connection between the Ito integral and Wick product 
is discussed, for example, in J^l Section 2.5]. 
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As an example, consider the following equation: 

(3.8) ut(t, x) = a(x)u xx (t, x) + b(x)u x (t, x) + u x (t, x) o W(t, x), < t < T, 
with initial condition u(0,x) = uq(x). In (|3.8[) . 

(WN1) W is the white noise process on R 2 . 

(WN2) The initial condition uq and the coefficients a, b are bounded and have continuous 

bounded derivatives up to second order. 
(WN3) There exists a positive number e so that a(x) > e, x G R. 
( WN4) The second-order derivative of a is uniformly Holder continuous. 

The equivalent Ito formulation of (|3.8j) is 

(3.9) du(t, x) = (a(x)u xx (t, x) + b(x)u x (t, x))dt + ek(x)u x (t, x)dwk(x), 
where {e/-, k > 1} is the Hermite basis in L2CR). 

With AikV = ei~v x , we see that condition 1)2. 2 Jl does not hold in any Sobolev space H^M). 
In fact, no traditional solution exists in any normal triple of Sobolev space. On the other 
hand, with a suitable definition of solution, equation 1)3.8(1 is solvable in the space (<S)_o of 
Hida distributions. 

Definition 3.10. A mapping u : R rf — > (S)- p is called weakly differentiable with respect to 
Xi at a point x* G R^ if and only if there exists a Ui{x*) G (S)- p so that, for all tp G (S) p , 
Di(u(x),ip)\x=x* = (Ui(x*),cp) . In that case, we write Ui{x*) = Diu(x*). 

Definition 3.11. A mapping u from [0, T] x R to (<S)_o is called a white noise solution of 
\S. 8\) if and only if 

(1) The weak derivatives ut, u x , and u xx exist, in the sense of Definition ^. H\ for all 
(t,x) G (0,T) x R. 

(2) Equality holds for all (t,x) G (0,T) x R d . 

(3) limjjo u{t, x) = uq(x) in the topology o/(<S)_o- 

Theorem 3.12. Under assumptions (WN1)-(WN4), there exists a white noise solution of 
US. 8\) . This solution is unique in the class of weakly measurable mappings v from (0,T) x R 
to (S)-o, for which there exists a non-negative integer q and a positive number K so that 




v(t, x)||_o,-? e dxdt < 00. 



Proof. Consider the S-transformed equation 

(3.10) F t (t,x;h) = a(x)F xx (t,x;h) + b(x)F x (t,x;h) + F x (t, x; h)h, 

< t < T, xGR, h E 5(R), with initial condition .F(0, a;; h) = uq(x). This a deterministic 
parabolic equation, and one can show, using the probabilistic representation of F, that 
F,Ft,F x , and F xx belong to U°. Then the inverse S-transform of F is a solution of (|3.8|) . 
and the uniqueness follows from the uniqueness for equation (|3.1U|) . The details of the proof 
are in jlU], where a similar equation is considered for x G R rf . □ 



Even though the initial condition in (|3.8|) is deterministic, there are no measurability re- 
strictions on no for the white noise solution to exist; see |12j for more details. 
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With appropriate modifications, the white noise solution can be defined for equations more 
general than l|3.8|) . The solution F = F(t, x; h) of the corresponding S-transformed equation 
determines the regularity of the white noise solution j!2l Section 4.1]. 

Two main advantages of the white noise approach over the Hilbert space approach are 

(1) no need for parabolicity condition; 

(2) no measurability restrictions on the input data. 

Still, there are substantial limitations: 

(1) There seems to be little or no connection between the white noise solution and the 
traditional solution. While white noise solution can, in principle, be constructed for 
equation (|2.7|) . this solution will be very different from the traditional solution. 

(2) There are no clear ways of computing the solution numerically, even with available 
representations of the Feynmann-Kac type |121 Chapter 4]. 

(3) The white noise solution, being constructed on a special white noise probability 
space, is weak in the probabilistic sense. Path-wise uniqueness does not apply to 
such solutions because of the "averaging" nature of the solution spaces. 



The objective of this section is to introduce the space of generalized random elements on 
an arbitrary stochastic basis. 

Let F = (Q, IF, {Ft}t>0: P) be a stochastic basis with the usual assumptions and Y, a 
separable Hilbert space with inner product (v)y and an orthonormal basis {y^, k > 1}. 
On F and Y, consider a cylindrical Brownian motion W, that is, a family of continuous 
^-adapted Gaussian martingales W y (t), y £ Y, so that W y (0) = and K(W yi (t)W y2 (s)) = 
min(i, s)(yi, 2/2) Y • In particular, 



are independent standard Wiener processes on F. 

Equivalently, instead of the process W, the starting point can be a system of independent 
standard Wiener processes {wk, k > 1} on F. Then, given a separable Hilbert space Y 
with an orthonormal basis {yk, k > 1}, the corresponding cylindrical Brownian motion W 
is defined by 



Fix a non-random T £ (0, 00) and denote by Jj! the sigma-algebra generated by Wk(t), k > 
1, < t < T. Denote by L2(W) the collection of Fjf -measurable square integrable random 
variables. 

We now review construction of the Cameron-Martin basis in the Hilbert space Z/2(W). 

Let m = {mk, k > 1} be an orthonormal basis in ^((0, T)) so that each belongs to 
Axj((0,T)). Define the independent standard Gaussian random variables 



4. Generalized Functions on the Wiener Chaos Space 



(4.1) 



W k (t)=Wy k (t), k>\, t>0, 



(4.2) 
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Consider the collection of multi-indices 

J={a = (a$, i,k>l), af £ {0, 1,2,.. .}, £}o*<oo}. 
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i . k 



The set J is countable, and, for every a £ J , only finitely many of a\ are not equal to zero. 
The upper and lower indices in oft represent, respectively, the space and time components 
of the noise process W. For a 6 J , define 



a 



Va! x . 



and 
(4.3) 



where H n is n th Hermite polynomial. For example, if 

/ 1 3 
2 4 
a = 000000 



with four non-zero entries a\ 



6,i 



4, then 

flsC^l) ^2(6,2) ^4fe,2) 



1; a 4 = 3; = 2; a 5 



V3! V2! V4! 
There are two main differences between (|3.2|) and (|4.3j) : 

(1) The basis Q4.3|) is constructed on an arbitrary probability space. 

(2) In (|4.3[) . there is a clear separation of the time and space components of the noise, 
and explicit presence of the time-dependent functions mi facilitates the analysis of 
evolution equations. 

Definition 4.1. The space ^(W) is called the Wiener Chaos space. The N-th Wiener 
Chaos is the linear subspace of L%(W), generated by £ a , \a\ = N. 

The following is another version of the classical results of Cameron and Martin [3]. 
Theorem 4.2. The collection H = a £ J} is an orthonormal basis in ^(W). 

We refer to S as the Cameron-Martin basis in L2(W). By Theorem 14. 21 every element v of 
L2(W) can be written as 

where v a = E(v£ a ). 

We now define the space £>(L2(W)) of test functions and the space T>'(L,2(W);X) of X- 
valued generalized random elements. 

Definition 4.3. 

(1) The space P(L2(W)) is the collection of elements from L2CW) that can be written in 
the form 
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for some v a £ R and a finite subset J v of J . 

(2) A sequence v n converges to v in 2?(L2(W)) if and only if J Vn C J v for all n and 
lim \v n a — v a \ = for all a. 

n—*oo ' 

Definition 4.4. For a linear topological space X define the space 

V(L2(W); X) of X -valued generalized random elements as the collection of continuous lin- 
ear maps from the linear topological space V(L2(W)) to X. Similarly, the elements of 
V (L2(W); Li((0,T); X)) are called X-valued generalized random processes. 

The element u of T>'(L2(W);X) can be identified with a formal Fourier series 

where u a £ X are the generalized Fourier coefficients of u. For such a series and for 
v £ V(L 2 (W)), we have 

u(v) = ^2 v » u »- 

Conversely, for u £ P'(L2(W); A), we define the formal Fourier series of u by setting 
u a = u(£ a ). If u £ L 2 (W), then u £ £>'(L 2 (W);R) and u(v) = E(uv). 

By Definition 14, 4( a sequence {u n , n > 1} converges to u in P'(L2(W); A) if and only if 
u n (v) converges to u(v) in the topology of X for every v £ T>(W). In terms of generalized 
Fourier coefficients, this is equivalent to lim u na = u a in the topology of X for every 

a £ J. 

The construction of the space D'(L2(W); X) can be extended to Hilbert spaces other than 
L2(W). Let H be a real separable Hilbert space with an orthonormal basis {e^, k > 1}. 
Define the space 

T)(H) = £ H : v = VkCk, Vk £ R, Jv — a finite subset of {1, 2, . . .} j . 

keJv 

By definition, v n converges to v in T>(H) as n — > 00 if and only if J Vn C J v for all n and 
lim \v n k — Vk\ =0 for all k. 

n— >oo ' 

For a linear topological space X, T>'(H; X) is the space of continuous linear maps from T)(H) 
to X. An element g of V'(H; X) can be identified with a formal series X]fc>i 3k ® Cfc so that 
5fc = 5(ejfc) G X and, for u G g(v) = J2kej v 9kVk- If X = R and J2k>idl < °o, then 

5 = ^2k>i9k e k £ and = (g,v)u, the inner product in ff. The space X is naturally 
imbedded into £>'(#; X): if u £ X, then £ fc>1 u <8> e k £ V'(H; X). 

A sequence g n = ^2 k >i 9n,k <8> e fc , n > 1, converges to 5 = Sfc>i 5fc ® e fc in A) if and 

only if, for every k > 1, lim g n k = 9k in the topology of A. 

n— »oo ' 

A collection {£fc, fc > 1} of linear operators from Ai to A2 naturally defines a linear 
operator C from V'(H;Xi) to V'(H;X 2 ): 

\fc>i / fc>i 

Similarly, a linear operator £ : T>'(H; X±) — > T>'(H; X2) can be identified with a collection 
A; > 1} of linear operators from Ai to A2 by setting Ck{u) = C(u (8> e&). Introduction 
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of spaces V'(H;X) and the corresponding operators makes it possible to avoid conditions 
of the type lt2~2l . 



5. The Malliavin Derivative and its Adjoint 

In this section, we define an analog of the ltd stochastic integral for generalized random 
processes. 

All notations from the previous section will remain in force. In particular, Y is a separable 
Hilbert space with a fixed orthonormal basis {y k , k > 1}, and S = a G J^}, the 

Cameron-Martin basis in L 2 (W) defined in (J4.3|) . 

We start with a brief review of the Malliavin calculus |37| . 

The Malliavin derivative B is a continuous linear operator from 

(5.1) 4(W) = |ii G L 2 (W) : H n " < °°} 

to L 2 (W; (L 2 ((0, T)) X Y)). In particular, 



(5-2) (B£ a )(t) = £ ^aUa-{i,k)mi(t)y, 



i.k 



where a (i, /c) is the multi-index with the components 



a (i, k) 



1 C max(a- —1,0), if i = j and k = I, 
■j L a j' otherwise. 



Note that, for each t G [0, T], B£ a (i) G £>(L 2 (W) x F). Using (f5"3|l . we extend the operator 
B by linearity to the space £>'(L 2 (W)): 



^ U a £a = Yl \ Ua Yl 
\a&J ) a<=J \ i.k 



For the sake of completeness and to justify further definitions, let us establish connection 
between the Malliavin derivative and the stochastic Ito integral. 

If u is an J^-adapted process from L 2 (W; L 2 ((0, T); Y)), then u(t) = Yl k >\ u k(t)yk, where 
the random variable u k (t) is J-j^-measurable for each t and k, and 



l-T 

V / E\u k (t)\ 2 dt < oo. 

rrf Jo 



k>l 

We define the stochastic Ito integral 

(5.3) U(t) = [\u( S ),dW(s)) Y = V f u k (s)dw k (s) 

k>i Jo 



o 



Note that U(t) is -measurable and E|[/(t)| 2 = ^ fc > x / Q * E\u k (s)\ 2 ds. 

The next result establishes a connection between the Malliavin derivative and the stochastic 
Ito integral. 
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Lemma 5.1. Suppose that u is an J-^ -adapted process from 

1/2 (W; 1/2 ((0, T); Y)), and define the process U according to \5.3\) . Then, for every < t < 
T and a € J , 

(5.4) E(*7(i)£ Q ) = E ! (u(s), (B^ a )(s)) Y ds. 

Jo 

Proof. Define = E(£ a | .T^). It is known (see or Remark EHH below) that 

(5.5) d£ a (t) = y<4€a-(i,k)(t)mi(t)dwk(t). 

i,k 

Due to J-^-measurability of Ufc(i), we have 

(5.6) u Ka (t) = E(« fc (t)E(Ca|^)) = E(u k (t)Ut))- 

The definition of U implies dU(t) = Ylk>i u k(t)dwk(t), so that, by (|5.5|l . (|5.6|) . and the ltd 
formula, 



(5.7) U a (t) = E(U(t)i a ) = [ y^Jc4u Ka - {i ^(s)m t (s)ds. 

Jo i,k 

Together with (|5.2|) . the last equality implies (|5.4j) . Lemma l5~Tl is proved. □ 

Note that the coefficients u k>a of u £ L 2 (W; L 2 ((0,T); H)) belong to L 2 ((0,T)). We there- 
fore define Uk, a ,i = f$ Uk,a{t)mi{t)dt. Then, by (|5.7j) . 



(5-8) l/ q (T) = X) V a i u k,<*-(i,k),i- 

Since U{T) = ^2 a& j U a (T)^ a , we shift the summation index in ()5.8|) and conclude that 



' *' ' " ' ( n (/./,) ) ^ ; otherwise. 



where 

+ 1, if z = j and k = I, 

As a result, t/(T) = 5(u), where 5 is the adjoint of the Malliavin derivative, also known as 
the Skorokhod integral; see jSZ] or [3S] f° r details. 

Lemma 15.11 suggests the following definition. For an //^-adapted process u from 
L 2 (W;L 2 ((0,T))), let W k u be the ^-adapted process from L 2 (W; L 2 ((0, T))) so that 

(5.11) (B* k u) a (t)=[ ^J^u a -a k )(s)mi(s)ds. 

Jo 

If u G L 2 (W; L2((0, T);Y)) is //"/^-adapted, then u is in the domain of the operator 5 and 
S{uI(s<t))=Z k >x(P>k)(t). 

We now extend the operators B£ to the generalized random processes. Let X be a Banach 
space with norm || • \\x- 
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Definition 5.2. If u is an X -valued generalized random process, then W k u is the X -valued 
generalized random process so that 

ft 



(5.12) (Bfew)a(i) = Y\ [ u a - (hk) {s)Jc^m l (s)ds. 

; JO 



If g G T>' \Y;T>' (L2(W); Lx((0,T); X)) j , then Wg is the X-valued generalized random 
process so that, for g = J2k>i9k ®Vk, 9k 6 T>'(L 2 (W); Li((0, T); X)), 

(5.13) (H*g) a (t) = V(D^fc) a (t) = V f 9 k ,a-{i,k){s)\f^ i rn i {s)ds. 

Using (j5.2j) . we get a generalization of equality (|5,4j) : 

(5.14) (B*«/) a (*)= [ g(BUs))(s)ds. 

Jo 

Indeed, by linearity, 

Theorem 5.3. IfT< oo, then B£ and B* are continuous linear operators. 



(J^mi( s )€a-{i,k)) ( s ) = \f4™i{s)g k:a - {iik) )(s). 



Proof. It is enough to show that, if u, u n G P' (Z^C? 7 ^); £i((0, T); X)) and 
linij^oo — u^qjII^^o.t)^) = f° r every a G J ', then, for every k > 1 and a 6 J, 
lim^oo ||(B*u) a - (B^ n ) a || Ll((0 ,T);X) = 0. 

Using (|5.12l) . we find 

\\(P%u) a ~ (D* k u n ) a \\ x (t) < V / \/af\\u a -, ik ) - u na -u k) \\x(s)\mi(s)\ds. 

i J ° 

Note that the sum contains finitely many terms. By assumption, | rrzj (t) | < Cj, and so 



%U) a ~ (n* k U n ) a \\ Ll ^ 0jT y t x) <C(a)^2y a^\\u a -^ - -W n ,a-(i,fc)llLi((0,T);X)- 

i 

Theorem 15.31 is proved. □ 



6. The Wiener Chaos Solution and the Propagator 



In this section we build on the ideas from |25| to introduce the Wiener Chaos solution and 
the corresponding propagator for a general stochastic evolution equation. The notations 
from Sections 0] and El will remain in force. It will be convenient to interpret the cylindrical 
Brownian motion W as a collection {w k , k > 1} of independent standard Wiener processes. 
As before, T G (0, oo) is fixed and non-random. Introduce the following objects: 

• The Banach spaces A, X, and U so that [/CI. 

• Linear operators 

A : Li((0,T);A) -» Li((0,T);X) and 
M k :L 1 ((0,T);A)^L 1 ((0,T);X). 

• Generalized random processes / G V (L2(W); Li((0, T); X)) and 
g k GP / (L 2 (W);L 1 ((0,T);X)). 
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• The initial condition no £ V (L2(W); U). 
Consider the deterministic equation 

(6.1) v (t)=v + [ {Av)(s)ds + [ ip{s)ds, 

Jo Jo 

where v £ U and ip £ Li((0, T);X). 

Definition 6.1. A function v is called a w(A,X) solution of i6.1\) if and only if v £ 
i>i((0, T);A) and equality holds in the space Li((0, T);A). 

Definition 6.2. An A-valued generalized random process u is called a w(A,X) Wiener 
Chaos solution of the stochastic differential equation 

(6.2) du{t) = (Au(t) + f(t))dt + (M k u(t) + g k (t))dw k {t), < t < T, u\ t=0 = u , 
if and only if the equality 

(6.3) u (t)=u + [ {Au + f)(s)ds + y^(Bl{M k u + g k ))(t) 

Jo k>i 

holds in V (L 2 (W);Li((0,T);X)). 

Sometimes, to stress the dependence of the Wiener Chaos solution on the terminal time T, 
the notation wt(A, X) will be used. 

Equalities (|6,3|) (|5.13|) mean that, for every a £ J , the generalized Fourier coefficient u a of 
u satisfies 

(6.4) u a (t) = n ,, 



+ j (Au + f) a (s)ds + [ Vi/of(M^ + 9k) a -(i,k)(s)mi(s)ds. 

J ° J ° Lk 



Definition 6.3. System \6.4\) is called the propagator for equation \6.£fy . 



The propagator is a lower triangular system. Indeed, If a = (0), that is, \a\ = 0, then the 
corresponding equation in (|6,4|) becomes 

(6.5) u (0) {t) = n ,( ) + I {Au (0) {s) + f (0) (s))ds. 

Jo 

If a = (j£), that is, oq = 1 for some fixed j and t and af = for all other i,k > 1, then the 
corresponding equation in ()6.4|) becomes 

u (je)(t) = uo,(ji) + / (Au {je) (s) + f( je) (s))ds 

(6.6) Jo 

+ / {Mkur )(s) + gu ) (s))mj(s)ds. 
Jo 

Continuing in this way, we conclude that (|6.4|) can be solved by induction on \a\ as long as 
the corresponding deterministic equation l)6.1j) is solvable. The precise result is as follows. 

Theorem 6.4. If, for every vq £ U and £ Li((0, T); X), equation h6.1\) has a unique 
w(A, X) solution v(t) = V(t, vq, ip), then equation \6. ty) has a unique w(A, X) Wiener Chaos 
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solution so that 

U a (t) = V(t,U 0t a,fa) + ^2 \f^ V ( t ^' m iMkU a -(i, k )) 

(6.7) 



+ Yl V a i V ^' °> m i9k,a-(i,k))- 



Proof. Using the assumptions of the theorem and linearity, we conclude that (|6.7I) is the 
unique solution of (|6.4|) , □ 

To derive a more explicit formula for u a , we need some additional constructions. For every 
multi-index a with \a\ = n, define the characteristic set K a of a so that 

K a = {(if , kf ), . . . , (i", &")}, 

ij < if < . . . < C and if if = then kf < kf +1 . The first pair (if, fcf ) in if a is the 
position numbers of the first nonzero element of a. The second pair is the same as the 
first if the first nonzero element of a is greater than one; otherwise, the second pair is the 
position numbers of the second nonzero element of a and so on. As a result, if a k > 0, then 
exactly pairs in K a are equal to (i, k). For example, if 

/0102300---\ 
1 2 1 ••• 

a = ooooooo--- 



V : : : : : : : 

with nonzero elements 

a\ = a\ = a\ = 1, a\ = a\ = 2, a\ = 3, 
then the characteristic set is 

K a = {(l,2), (2,1), (2,2), (2,2), (4,1), (4,1), (5,1), (5,1), (5,1), (6,2)}. 
Theorem 6.5. Assume that 

(1) for every vq £ U and if £ L\((0, T);X), equation \ti '. 1\) has a unique w(A, X) 
solution v(t) = V(t, vo,ip), 

(2) the input data in j6.4\ ) satisfy = and f a = uo, a = if \a\ > 0. 

Let U(o)(i) = V(t,uo,0) be the solution of \6.4\) for \a\ = 0. For a £ J with \a\ = n > 1 
and the characteristic set K a , define functions F n = F n (t; a) by induction as follows: 

F l {t-a) = V(t,0, mi M k u m ) if K a = {(i,k)}; 

n 

(6.8) F n (t;a) = ^2v(t,0,m ij M kj F n - 1 (-,a~(i J ,k,))) 

i=i 

if K a = {(ix,kx),...,{i n ,k n )}. 

Then 

(6.9) u a {t) = ^=F n {t-a). 

ya\ 



Proof. If \a\ = 1, then representation (J6.9|) follows from (|6.6|) . For |a| > 1, observe that 
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• If u a (t) = \fa\u a and \a\ > 1, then (|6.4|) implies 

u(t) = / Au a {s)ds + E / a.imi(s)MkU a -(i,k)(s)ds. 
Jo i,k Jo 

• If K a = {(it, k%), ... ,(i n , k n )}, then, for every j = l,...,n, the characteristic set 

-(ij,kj) °^ a fe> %) * s obtained from by removing the pair (ij, kj). 

• By the definition of the characteristic set, 

n 

^aimi(s)M k u a - {itk) (s) = ^2m ij (s)M kj u a - {ij ^ : . ) (s). 
i,k j=l 

As a result, representation (|6.9() follows by induction on \a\ using (|6.7|) : 
if |a| = n > 1, then 

n 

= E^'°' m V^ fe ^ a ~fe- fc i)) 

(6.10) 



= E V & °» ^Xfe-F^-^C-; Q-(i j5 %)) = a). 
i=i 

Theorem 16.51 is proved. □ 

Corollary 6.6. Assume that the operator A is a generator of a strongly continuous semi- 
group <I> = <&t,s, t > s > 0, in some Hilbert space H so that A C H , each M k is a bounded 
operator from A to H, and the solution V(t, 0, ip) of equation hb\l\) is written as 

(6.11) V(t,0,<p)= f $tMs)ds, <peL 2 ({0,T);H)). 

Jo 

Denote by V n the permutation group of {1, . . . , n}. Ifur \ £ ^((0, T); H)), then, for \a\ = 
n > 1 with the characteristic set K a = {(ii, k\), . . . , (i n , kn)}, representation \6. y\) becomes 



(6.12) 



1 , /"* f Sn f S2 
VcA ^ n Jo Jo Jo 

<5>t,s n M ka{n) ■ ■ ■ $ S2 , Sl A4fc CT(1) u (0 )(si)m V(n) (s n ) • • • m ia{1) (*i)dsi . . . ds n . 



Also, 

ft f-S n f-S2 



(6-13) \ a \ =n kl ,...,k n >i Jo Jo Jo 

$t,s n Mk n ■••* S2)S1 (M kl u^ + g kl (s 1 )) dw kl (si) ■ ■ ■ dw kn (s n ), n > 1, 

and, for every Hilbert space X, the following energy equality holds: 

_°° rt rs n r-s2 

Eim*)Hx= E //■■■/ 

(6-14) k u lX n =l J ° J ° J ° 

\\^t,s„M kn • • • ^s^M^u^is^Wxdsx . . . ds n ; 
both sides in the last equality can be infinite. For n = 1, formulas \6.1ty) and become 



(6.15) = / ®t,sMkU( )(s) mi(s)ds; 
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oo . t 

(6.16) X) ll««(*)llx=Z)/ \\®t,sM k u i0 )(s)\\ x ds. 

\a\=l k=l J ° 

Proof. Using the semi-group representation (|6.11|) . we conclude that 1)6.12(1 is just an ex- 
panded version of (|6.9|) . 

Since {rrii, i > 1} is an orthonormal basis in £2(0, T), equality (|6.16|) follows from (|6.15|) and 
the Parcevall identity. Similarly, equality (|6.14() will follow from (|6.12|) after an application 
of an appropriate Parcevall's identity. 

To carry out the necessary arguments when \a\ > 1, denote by J\ the collection of one- 
dimensional multi-indices (3 = /?2, ■ ■ ■) so that each is a non-negative integer and 
\(3\ = Y^,i>i A < 00 • Given a /3 G Ji with |/?| = n, we define Kp = {ii,...,i n }, the 
characteristic set of (3 and the function 

(6-17) E p (si,...,s n ) = ^ n»ii(sa(i)) --,m <n(M»))- 

By construction, the collection {Ep,f3 £ J\, |/?| = n} is an orthonormal basis in the sub- 
space of symmetric functions in I/2((0, T) n ; X). 

Next, we re-write (|6.12() in a symmetrized form. To make the notations shorter, denote 
by s( n ) the ordered set (si,...,s n ) and write ds n = ds\...ds n . Fix t £ (0,T] and the 
set fc( n ) = {ki, . . . , k n } of the second components of the characteristic set K a . Define the 
symmetric function 

G(t,A; (n) ;s (n) ) 

lU ' l " i = ~T= } Y^ ^.^(n)- M ^---^s CT(2) , S(T(1) ^fc 1 U ( o)(s (J (l))ls (T(1) <...< S(T(n) <t(s (n) ). 

vn. aeVn 
Then (|6.12|) becomes 

(6.19) ««(*)=/ G(t,fc(">; S < n >)^ (a) (>>)ds» 

where the multi-indices a and /3(a) are related via their characteristic sets: if 

K a = {(ii,ki), . . . , (i n , k n )}, 

then 

K 0(a) = {h, ■ ■ ■ ,in}- 

Equality (|6.19|) means that, for fixed k^ n \ the function u a is a Fourier coefficient of the 
symmetric function G(t, k^; s^) in the space L2((0,T) n ; X). Parcevall's identity and 
summation over all possible k^ yield 

which, due to (|6.18|) . is the same as ()6.14|) . 

To prove equality (|6.13|) . relating the Cameron-Martin and multiple ltd integral expansions 
of the solution, we use the following result |13l Theorem 3.1]: 

1 l' Sn l' S2 

Ca = —F^ / ••• / ^/3(a)(s (n) )d w fci( s l) • • -dWk^Sn); 



'a\ Jo Jo 
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see also (3Z1 PP- 12-13]. Since the collection of all Ep is an orthonormal basis, equality 
(|6.13j) follows from (|6.19j) after summation over al k\ , . . . , kn . 

Corollary 16.61 is proved. □ 

We now present several examples to illustrate the general results. 
Example 6.7. Consider the following equation: 

(6.20) du(t, x) = (au xx (t, x) + f(t, x))dt + (au x (t, x) + g(t, x))dw(t), t > 0, x £R, 

where a > 0, a G R, / G L 2 ((0, T); R 2 X (R)), o £ L 2 ( (0, T); L 2 (M)), and u| t=0 = «o 6 
L 2 (M). By Theorem E31 if cr 2 < 2a, then equation (|6~2TH) has a unique traditional solution 
uGL 2 (W;L 2 ((0,T); J ff 2 1 (lR))). 

By J^-measurability of u(t), we have 

E(«(t)£ a ) = E(u(t)E(&|^)). 
Using the relation (|5,5|) and the Ito formula, we find that u a satisfy 

du a = a(u a ) xx dt + ^ \fa~icr(ua-(i))xmi(t)dt, 

i 

which is precisely the propagator for equation (|6.2UI) . In other words, if 2a > a 2 , then the 
traditional solution of (|6.2U|) coincides with the Wiener Chaos solution. 

On the other hand, the heat equation 

v (t, x) = vq(x) + / v xx (s,x)ds + / tp(s,x)ds, vo G -L 2 (M) 
Jo Jo 

with ip G L 2 ((0, T); if^M)) has a unique w(H%(R), H^iR)) solution. Therefore, by Theo- 
rem !6.41 the unique ui(-£f 2 (IR), fl^" (M)) Wiener Chaos solution of ()6.2(Jj) exists for all cr G R. 

In the next example, the equation, although not parabolic, can be solved explicitly. 
Example 6.8. Consider the following equation: 

(6.21) du(t,x) = u x (t,x)dw(t), t > 0, x G R; u(0,x)=x. 
Clearly, u{t,x) = x + w(t) satisfies (|6,21[) , 

To find the Wiener Chaos solution of ()6.21j) . note that, with one-dimensional Wiener process, 
a\ = on, and the propagator in this case becomes 

u a (t,x) = xl(\a\ = 0) + / y/ai(u a -(i)(s,x)) x mi(s)ds. 

Jo ■ 

Then u a = if \a\ > 1, and 

(6.22) u(t,x) = x + / rrii(s)ds = x + w{t). 

i>i 

Even though Theorem E31 does not apply, the above arguments show that u(t, x) = x + w(t) 
is the unique w(A, X) Wiener Chaos solution of (|6.21j) for suitable spaces A and X, for 
example, 

X = if:J(l + x 2 )~ 2 f 2 (x)dx < oo | and A = {/:/,/' G X}. 
Section El provides a more detailed analysis of equation (|6.2ip . 
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If equation (|6.2|) is anticipating, that is, the initial condition is not deterministic and/or 
the free terms /, g are not J-^-adapted, then the Wiener Chaos solution generalizes the 
Skorohod integral interpretation of the equation. 

Example 6.9. Consider the equation 

(6.23) du(t, x) = -u xx (t, x)dt + u x (t, x)dw(t), t G (0, T], x G R, 
with initial condition u(0, x) = x 2 w(T). Since w(T) = VT^i, we find 

(6.24) (u a ) t (t,x) = -(u a ) xx (t,x) + Y^V^i m i( t )( u a-(i))x(t,x) 

i 

with initial condition u a (0, x) = \^Tx 2 I(\a\ = l,at\ = 1). By Theorem 16.41 there exists a 
unique w(A, X) Wiener Chaos solution of (|6.23|) for suitable spaces A and X. For example, 
we can take 

X = {/ : J (1 + x 2 )- 8 f(x)dx < ooj and A = {/ : /, /', /" G X}. 

System (j6.24[) can be solved explicitly. Indeed, u a = if \a\ = or \a\ > 3 or if a± = 0. 
Otherwise, writing Mj(t) = f~ rrii(s)ds, we find: 

u a (t, x) = (t + x 2 )VT, if \a\ = 1, a± = 1; 

u a (t, x) = 2V2 xt, if \a\ =2, «i = 2; 

u a (t, x) = 2\fT xMi(t), if \a\ =2, «i = a.{ = 1, 1 < i; 

u a (t,x) = y — t 2 , if |a| =3, a± = 3; 



u a {t,x) = 2V2T Mi(i)Mi(i), if |a| =3, ai = 2, a* = 1, 1< i; 

u a (t,x) = V2T M 2 (t), if |a| =3, a x = 1, a» = 2, 1 < i; 

u a (i,x) = 2-v/T Mi(t)Mj(t), if |a| =3, ai = a» = a 3 - = 1, 1 < i < j. 

Then 

(6.25) u(t, x) = = w ( r )^ 2 (0 - 2tw(t) + 2(W(T)«;(t) - t)x + x 2 w(T) 

is the Wiener Chaos solution of (|6.23|) . It can be verified using the properties of the Skorohod 
integral [HZJ that the function u defined by (|6.25j) satisfies 

u(t, x) = x 2 w(T) + - I u xx (s,x)ds+ I u x (s,x)dw(s), t&[0, T], 
2 Jo Jo 

where the stochastic integral is in the sense of Skorohod. 



7. Weighted Wiener Chaos Spaces and S-Transform 



The space P'(L2(W); X) is too big to provide any reasonable information about regularity of 
the Wiener Chaos solution. Introduction of weighted Wiener chaos spaces makes it possible 
to resolve this difficulty. 

As before, let H = a G J} be the Cameron-Martin basis in L2(W), and V(L2(W); X), 
the collection of finite linear combinations of £ a with coefficients in a Banach space X. 
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Definition 7.1. Given a collection {r a , a E J} of positive numbers, the space 1ZL2(W; X) 
is the closure of T>(L2(W); X) with respect to the norm 



2 _ 2 II ||2 

TZL 2 (W;X) ■— 2-^i r a.\\ v o.\\X- 

aej 



The operator 1Z defined by {lZv) a := r a v a is a linear homeomorphism from 1ZL2(W; X) to 
L 2 (W;X). 

There are several special choices of the weight sequence 1Z = {r a , a E J} and special 
notations for the corresponding weighted Wiener chaos spaces. 

• If Q = {qi, q2, . . .} is a sequence of positive numbers, define 

i,k 

The operator 1Z, corresponding to r a = q a , is denotes by Q. The space QL2(W; X) 
is denoted by L2,q(W; X) and is called a Q-weighted Wiener chaos space. The 
significance of this choice of weights will be explained shortly (see, in particular, 
Proposition I7,4|) . 

• If 

^ = («o p II( 2 ^) 7aifc ' p^ £R > 

i,k 

then the corresponding space TZL2(W; X) is denoted by (5) Pi7 (X). As always, the 
argument X will be omitted if X = R. Note the analogy with Definition 13,31 

The structure of weights in the spaces L2,q and (5) P , 7 is different, and in general these 
two classes of spaces are not related. There exist generalized random elements that belong 
to some L2,q(W;X), but do not belong to any (S) pa (X). For example, u = Ylk>i ek 

belongs to L 2! q(W) with q^ = e~ 2fc2 , but to no («S) P)7 , because the sum Y,k>i e 2k2 (k\) p (2k) 1 
diverges for every p, 7 G R. Similarly, there exist generalized random elements that belong 
to some (S)p^(X), but to no L 2i q(W; X). For example, u = Y^ n >i V^C(n)i where (n) is the 
multi-index with a\ = n and oq = elsewhere, belongs to (S)—\ _i, but does not belong to 
any L2,q(W), because the sum Yln>i diverges for every q > 0. 

The next result is the space-time analog of Proposition 2.3.3 in \1'2\ . 

Proposition 7.2. The sum 

a&J i,fc>l 

converges if and only if 7 > 1 . 
Proof. Note that 

(7.1) £ n = n ( E(( 2 ^)- 7 ) n ) = n a-L^y ? > ° 

aeji,k>l i,k>l \n>0 J i,k V y ' ' 

The infinite product on the right of (|7.1|) converges if and only if each of the sums Yli>i * _7 > 
Sfc>i ^ 7 converges, that is, if an only if 7 > 1. □ 

Corollary 7.3. For every u E T>'(W;X), there exists an operator 7Z so that 
TZu E L 2 (W;X). 
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Proof. Define 

.2 _. 1 TT fr,-u\-2a? 



Then 

'\u 



L + \\ u *\\x itK>1 



^iii 2( w ; x) = e itSv n w 20 * * e n < ~- 

oej " a|lj5<r i,fe>l aej i,k>l 

□ 



The importance of the operator Q in the study of stochastic equations is due to the fact 
that the operator 1Z maps a Wiener Chaos solution to a Wiener Chaos solution if and only 
TZ = Q for some sequence Q. Indeed, direct calculations show that the functions u a ,a £ J , 
satisfy the propagator (|6,4|) if and only if v a = (JZu) a satisfy 

v a (t) = (Ku ) a + I {Av + Kf) a (s)ds 
Jo 

(7-2) f t j- pa 

Therefore, the operator 1Z preserves the structure of the propagator if and only if 

Pa 

= Qk, 

Pa~(i,k) 

that is, p a = q a for some sequence Q. 

Below is the summary of the main properties of the operator Q. 
Proposition 7.4. 



(1) If qk < q < 1 for all k > 1, then L2,q(W) C (<S)o,- 7 for some 7 > 0. 

(2) If qk > q > 1 for all k, then L2,q(W) C L 2 l (W) for all n>l, that is, the elements 
of Ii2 t Q (W) are infinitely differentiate in the Malliavin sense. 

(3) If u 6 L2,q(W;X) with generalized Fourier coefficients u a satisfying the propaga- 
tor \6.4\ ), and v = Qu, then the corresponding system for the generalized Fourier 
coefficients of v is 

V a {t) = (Qu ) a + [ (Av + Qf) a {s)ds 
Jo 

(7.3) ft 

+ / ^y a i( M kV + Qgk) a -(i,k)(s)qkmi(s)ds. 

(4) The function u is a Wiener Chaos solution of 



(7.4) u(t) = u + / (Au(s) + f(s))dt + / (Mu(s) + g(s),dW(s)) Y 

Jo Jo 

if and only if v = Qu is a Wiener Chaos solution of 

(7.5) v(t) = (Qu) + [\av(s) + Qf(s))dt + f (Mv{s) + Qg(s), dW Q (s)) Y , 

J Jo 

where, for heY, W®(t) = Ek>l(h ,yk)Yqkw k {t). 
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The following examples demonstrate how the operator Q helps with the analysis of various 
stochastic evolution equations. 

Example 7.5. Consider the w(H 2 (R), -ff^~ 1 (M)) Wiener Chaos solution u of equation 

(7.6) du(t, x) = (au xx (t, x) + f(t, x))dt + au x (t, x)dw(t), < t < T, i£R, 

with / G L 2 (tt x (0,T);H 2 1 (R)), g G L 2 (Q x (0, T); L 2 (R)), and u|t=o = u £ L 2 (R). 
Assume that a > and define the sequence Q so that qk = q for all > 1 and q < v2a/er. 
By Theorem 12 .5| equation 

dv = (av xx + f)dt + (qou x + o)dw; 

with v\t=o = uq, has a unique traditional solution 

v e L 2 (W;L 2 ((0,T); J ff 2 1 (M))) f]L 2 (W; C((0, T); L 2 (R))) . 

By Proposition 17.41 the w{H\ (R), ff 2 ~ 1 (R)) Wiener Chaos solution u of equation (|7.6|) sat- 
isfies u = Q _1 f and 

« G ^2,Q (W; L 2 ((0, T); ^(R))) f| L 2 , Q (W; C((0, T); L 2 (R))) . 

Note that if equation Q7.6|) is strongly parabolic, that is, 2a > c 2 , then the weight g can be 
taken bigger than one, and, according to the first statement of Proposition 17.41 regularity 
of the solution is better than the one guaranteed by Theorem 12. 51 

Example 7.6. The Wiener Chaos solutions can be constructed for stochastic ordinary 
differential equations. Consider, for example, 

(7.7) u{t) = 1+ / J^u{s)dw k {s), 

Jo k>i 

which clearly does not have a traditional solution. On the other hand, the unique to(R, R) 
Wiener Chaos solution of this equation belongs to L 2j q (W; L 2 ((Q, T)) for every Q satisfying 
J2k1k < 00 • Indeed, for (|7.7j) . equation (|7.5f) becomes 

v(t) = l+ / y2v(s)q k dw k (s). 
Jo k 

If ^2 k qf, < co, then the traditional solution of this equation exists and belongs to 
L 2 (W;L 2 ((0,T))). 

There exist equations for which the Wiener Chaos solution does not belong to any weighted 
Wiener chaos space L 2) q. An example is given below in Section ITU 

To define the S-transform, consider the following analog of the stochastic exponential HM.fijl . 
Lemma 7.7. If h G V (L 2 ((0, T); Y)) and 

£{h) = exp n\h(t),dW(t))y \\h(t)f Y dt) , 

then 

• £(h) G L 2> q(W) for every sequence Q. 

• £{h) G (S) pa for < p < 1 and 7 > 0. 

• £ (h) G (<5)i )7 , 7 > 0, as long as ||^||x 2 ((o t)-y) ^ s sufficiently small. 
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Proof. Recall that, if h G V(L,2((0,T);Y)), then h(t) = Yli kei h hk^imiifyVki where Ih is a 
finite set. Direct computations show that 



i,k \n>0 ' / aej V 



where h a =\\ ik h^. In particular, 
(7.8) (£(h)) a ''" 



Consequently, for every sequence Q of positive numbers, 
(7-9) \\£(h)\\l 2Q{w) = exp j hid ) < do. 

Similarly, for < p < 1 and 7 > 0, 

(7.10) iigwnk„ = En J,ni- = n e 



QSJ i,k v * 7 j,fce/ h \n>0 v 7 



< 00, 



and, for p = 1, 



(7.11) wmw's)^ = E II(( 2 ^) 7 M 2a ' = II E« 2 * fc ) 7 M 2n < °°> 

a&J i,k i,kelh \n>0 J 

if 2 (max( m n ) g / h )(mn) 7 ) £^ i < 1. Lemma [?~71 is proved. □ 

Remark 7.8. It is well-known (see, for example, 24, Proof of Theorem 5.5]) that the 
family {£(h),h G V (L2((0,T);Y))} is dense in L2(W) and consequently in every L 2j q(W) 
and every (5) Pi7 , — 1 < p < 1, 7 G M. 

Definition 7.9. If u € L 2 ,q{W;X) for some Q, or if u G U g >o( 5 )-P,-7pO> < P < 1 > 
then the deterministic function 

(7.12) Su(h) = £ ^ g x 

is called the S-transform of u. Similarly, for g G T>' (Y; L2 j q(W; X)) the S-transform 
Sg(h) G V'(Y;X) is defined by setting (Sg(h))k = (Sgk)(h). 



Note that if u G L 2 (W;X), then Su(/t) = E(u£(h)). If it belongs to L 2:Q (W;X) or to 
U g > (5)-p- 7 (X), < p < 1, then Su(h) is defined for all h G £> (L 2 ((0, T); Y)) . If n G 
|J 7>0 (5)_i i _ 7 (X), then Su(h) is defined only for /i sufficiently close to zero. 

By Remark l7.8l an element u from L2,q(W;X) or (J 7>0 (5)_ Pi _ 7 (X), < p < 1, is uniquely 
determined by the collection of deterministic functions Su(h), h G V (L2((0, T); Y)) . Since 
£ (/i) > for all h G 2? (L2((0, T); Y)), Remark 17.81 also suggests the following definition. 

Definition 7.10. An element u from L2,q(W) or U 7>0 (5)_ Pi _ 7 , < p < 1 is called non- 
negative (u>0) if and only if Su(h) > for all h G V (L2((0, T); 1")). 



The definition of the operator Q and Definition 17. 1UI imply the following result. 
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Proposition 7.11. A generalized random element u from L2,q(W) is non-negative if and 
only if Qu > 0. 

For example, the solution of equation (|7.7|) is non-negative because 

Qu(t) = exp ^2(q k w k (t) - (l/2)gg) 

\k>l 

We conclude this section with one technical remark. 

Definition 17.91 expresses the S-transform in terms of the generalized Fourier coefficients. 
The following results makes it possible to recover generalized Fourier coefficients from the 
corresponding S-transform. 

Proposition 7.12. If u belongs to some L2,q(W; X) or U 7>0 (5)_ P; _ 7 (A") ; < p < 1, then 

(T.i3) „„ =i ( n «W 



()l>'i,, 



h=0 



Proof. For each a G J with K non-zero entries, equality (|7.12|) and Lemma l7.7l imply that 
the function Su(h), as a function of K variables h k> i, is analytic in some neighborhood of 
zero. Then (|7.13j) follows after differentiation of the series (|7,12jl . □ 

8. General Properties of the Wiener Chaos Solutions 

Using notations and assumptions from Section H3 consider the linear evolution equation 

(8.1) du(t) = (Au(t) + f{t))dt + (Mu(t) + g(t), dW{t)) Y , < t < T, u\ t=0 = u . 

The objective of this section is to study how the Wiener Chaos compares with the traditional 
and white noise solutions. 

To make the presentation shorter, call an X-valued generalized random element S-admissible 
if and only if it belongs to L2 j q(J 7W ; X) for some Q or to (S) Pjq (X) for some p G [—1, 1] 
and gGR. It was shown in Sectional that, for every S-admissible u, the S-transform Su(h) 
is defined when h = Y^i k ^k-i^iVk £ ^(-^((O, T); Y)) and is an analytic function of h^i in 
some neighborhood of h = 0. 

The next result describes the S-transform of the Wiener Chaos solution. 
Theorem 8.1. Assume that 

(1) there exists a unique w(A,X) Wiener Chaos solution u of \8.1\) and u is S- 
admissible; 

(2) For each t G [0, T], the linear operators A(t), Aik(t) are bounded from A to X; 

(3) the generalized random elements uo,f,gk are S-admissible. 

Then, for every h G V(L2((0, T); Y)) with |H|^ 2 (( T)-Y) sufficiently small, the function 
v = Su(h) is a w(A, X) solution of the deterministic equation 

(8.2) v(t) = Su (h) + / (Av + Sf(h) + (M k v + Sg k (h))h k ^j (s)ds. 
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Proof. By assumption, Su(h) exists for suitable functions h. Then the S-transformed equa- 



tion (|8.2j) follows from the definition of the S-transform (|7.12j) and the propagator equation 
(|6.4j) satisfied by the generalized Fourier coefficients of u. Indeed, continuity of operator A 
implies 

h a , h° 



S{Au)(h) = S2 —j=Au a = AS^^=u a =A(Su(h)). 
~f Vd ~? va! 



Similarly, 



~T= x Yl V a> i M k u a-{i,k) m i = Y Y I -(■ ^7 M kU a -{i,k) m i h k 



= ^ | y~] -j=M k u a mih kti = M k {Su(h))h k . 
i,k V a * / 

Computations for the other terms are similar. Theorem 18. II is proved. 
Remark 8.2. If he V(L 2 ({0,T);Y)) and 



□ 



5.3) £ t {h) = exp (j 



(h(s),dW(s)) 



Y 



\f wmw 



ydt 



then, by the ltd formula, 

(8.4) d£ t {h)=£ t {h){h(t),dW{t)) Y . 

If uq is deterministic, f and g k are -adapted, and u is a square-integrable solution of 
18.1)) . then equality \8. 6 J\) is obtained by multiplying equations \8-4\l an d H8.1)) according to 
the ltd formula and taking the expectation. 

Remark 8.3. Rewriting \8.4\) as 

d£ t (h) = £t(h)h kti mi(t)dw k (t) 

and using the relations 

1 



£ t (h)=E(£ T (h)\^), e 




d a *£ T {h) 



dh 



k.i 



h=0 



we arrive at representation \5.5\) for E(^ Q ,|J r t w/ ). 



A partial converse of Theorem 18.11 is that, under some regularity conditions, the Wiener 
Chaos solution can be recovered from the solution of the S-transformed equation (|8.2|) . 

Theorem 8.4. Assume that the linear operators A(t), Ai k (t), t € [0, T], are bounded from 
A to X, the input data uq, f, g k are S-admissible, and, for every h e T>(L2((0, T); Y)) 
with ||^|lx 2 ((o t)-f) sufficiently small, there exists a w(A,X) solution v = v(t;h) of equation 
i8.ty) . We write h = h k ^miy k and consider v as a function of the variables h k ^. Assume 
that all the derivatives of v at the point h = exists, and, for a S J , define 



1.5) 



u a (t) 




d a *v(t;h) 



dh 



k:i 



h=0 



Then the generalized random process u(t) = Y2aeJ Ua ^^ a * s a W {A,X) Wiener Chaos 
solution of 18.1)) . 
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Proof. Differentiation of (|8.2|) and application of Proposition 17. 12l show that the functions 
u a satisfy the propag ator ([Ojl . □ 

Remark 8.5. The central part in the construction of the white noise solution of \8.1]) 
is proving that the solution of \8.lfy is an S-transform of a suitable generalized random 
process. For many particular cases of equation \8.1]) . the corresponding analysis is carried 
out in |10| 1121 I33[ I4f)j . The consequence of Theorems \8.1\ and \8.4\ is that a white noise 
solution of \8.1\) . if exists, must coincide with the Wiener Chaos solution. 

The next theorem establishes the connection between the Wiener Chaos solution and the 
traditional solution. Recall that the traditional, or square-integrable, solution of (|8.1|) was 
introduced in Definition 12.21 Accordingly, the notations from Section [2] will be used. 

Theorem 8.6. Let (V, H, V') be a normal triple of Hilbert spaces. Take a deterministic 
function uq and -adapted random processes function, f and g k so that \2. S\) holds. Under 
these assumptions we have the following two statements. 

(1) An J- -adapted traditional solution of \8.1]) is also a Wiener Chaos solution. 

(2) If u is a w(V, V) Wiener Chaos solution of \8.1\) so that 

(8.6) V f / \\u a (t)\\ 2 v dt+ sup <oo, 

then u is an -adapted traditional solution of \8.1]) . 

Proof. (1) If u = u(t) is an J^-adapted traditional solution, then 

u a (t) = E(u(t)Ca) = E {u(t)E(£ Q \ry )) = E(u(i)£«(t)). 

Then the propagator ()6.4|) for u a follows after applying the Ito formula to the product 
an d using (|5.5() . 

(2) Assumption ()8.6|) implies 

u E L 2 (n x (0, T); V) f| L 2 (ft; C((0, T); H)). 

Then, by Theorem 18.11 for every <p £ V and h £ T>((0,T);Y), the S-transform of u 
satisfies 

(uh(t), <p)h = {uo,<p)h + f (Au h (s),(p)ds + / (f(s),tp}ds 

Jo Jo 

+ ^"5^ / \/ a i m i( S ){( M k u a-(i,k){s),(p) H 
a£j °" i,k J ° 

+ (gk{s),cp) H I(\a\ = l))ds. 
If I(t) = f (Mku(s),<p)Hdw k (s), then 



(8.7) E(I(t)Ut))= [ Yj^ m ^ M k u a -{i,k){s)^) H ds. 

Jo i,k 

Similarly, 

E Ka(0 / (gk(s),(p) H dw k (s) ) =y~] / Ja>{mi(s)(g k (s),(p) H I(\a\ = i)ds. 
V Jo J - k Jo 
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Therefore, 



/i a , r t i — 

5Z — T / V a i m i( s ) {M. k u a - (i fc) (s) , ds 

al , JO 



aej ' i,k 

As a result, 

E (S(h)(u(t), <p) H ) = E (£(h)(u , <p)h) 

} + E f £(/») y ¥>>ds) + E ^(/i) y </(s), ^ds 

+ E f y ((Alfcu(a), <p)h + (9k(s), <p) H ) dw k (s) 

Equality (|8.8|) and Remark l7.8l implv that, for each t and each tp, (|2.4j) holds with probability 
one. Continuity of u implies that, for each ip, a single probability-one set can be chosen for 
all t £ [0, T]. Theorem 19.61 is proved. □ 

9. Regularity of the Wiener Chaos Solution 

Let F = J 7 , {^t}t>0jlP) be a stochastic basis with the usual assumptions and w k = 
w k (t), k > 1, t > 0, a collection of standard Wiener processes on F. As in Section [21 let 
(V,H,V) be a normal triple of Hilbert spaces and A(t) : V —* V, .Mfc(i) :V^>H, linear 
bounded operators; t £ [0, T]. 

In this section we study the linear equation 

(9.1) u (t) = u + [ (Au(s) + f{s))ds + [ {M k u{s) + g k (s))dw k , < t < T, 

Jo Jo 

under the following assumptions: 

Al There exist positive numbers C\ and 5 so that 

(9.2) (A(t)v,v) +S\\v\\l < CiWvfjj, v € V, t € [0,T]. 
A2 There exists a real number C2 so that 

(9.3) 2{A(t)v, v) + ^2 \\M k (t)vf H < C 2 \\v\\ 2 H , v € V, t € [0,T]. 

k>l 

A3 The initial condition uq is non-random and belongs to H; the process / = f(t) is 
deterministic and f Q \\f(t)\\y,dt < 00; each g k = g k {t) is a deterministic processes 

and E fc >i Jo llfffc(*)lllf^ < °°- 

Note that condition (|9.3f) is weaker than (|2.5|) . Traditional analysis of equation (|9.1|) under 
()9.3j) requires additional regularity assumptions on the input data and additional Hilbert 
space constructions beyond the normal triple 1421 Section 3.2]. In particular, no existence of 
a traditional solution is known under assumptions A1-A3, and the Wiener chaos approach 
provides new existence and regularity results for equation (|9.1|) . A different version of the 
following theorem is presented in |29| . 

Theorem 9.1. Under assumptions Al A3, for every T > 0, equation h9.1)) has a unique 
w(y,V) Wiener Chaos solution. This solution u = u(t) has the following properties: 
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(1) There exists a weight sequence Q so that 

u E L 2 , Q (W; L 2 ((0, T); V)) f] L 2 ,q(W; C((0, T); H)). 

(2) For every 0<t<T, u(t) E L 2 (0; H) and 

(9.4) E||u(t)|& < 3e c ^ ( \\u f H + C f f \\f{s)\$,ds + £ f bk{s)f H ds 

Jo k>1 Jo 



where the number C 2 is from \9.3\) and the positive number Cf depends only on 5 
and C\ from \9. #|) . 
(3) For every < t <T, 

ft rsn rs2 

<t) = u(o) + E E //•••/ 

(9-5) n>l kl ,...,k n >l J ° J ° J ° 

$t,s n M kn ■ ■ ■ $s 2 , Sl (M kl U( 0) +g kl (si)) dw kl (s!) ■ ■■dw kn (s n ), 
where <&t,s is the semi-group of the operator A. 

Proof. Assumption A2 and the properties of the normal triple imply that there exists a 
positive number C* so that 

(9.6) J2 WM&Mh < C*\\v\\ 2 v , veV, t€ [0,T]. 

fc>l 

Define the sequence Q so that 

(9-7) q k = (j^) 12 :=q, k > 1, 

where \i E (0, 2) and 5 is from Assumption Al. Then, by Assumption A2, 

(9.8) 2{Av,v) + £ q 2 \\M k v\\ 2 H < -(2 - v)5\\v\\ 2 v + Ci\\v\\ 2 H . 

k>\ 

It follows from Theorem 12.41 that equation 

(9.9) v{t) =u + [ (Av + f)(s)ds + Y] f q{M k v + g k )(s)dw k (s) 

Jo k ^Jo 

has a unique solution 

v E L 2 (W; L 2 ((0, T); V)) f| L 2 (W; C((0, T); H)). 

Comparison of the propagators for equations ()9.1|) and (|9.9|) shows that u = Q~ l v is the 
unique w(V,V) solution of (|9.1|) and 

(9.10) u E L 2 , Q (W; L 2 ((0, T); F)) f| L 2 , Q (W; C((0, T); #)). 



If C* < 25, then equation ()9.1j) is strongly parabolic and q > 1 is an admissible choice of 
the weight. As a result, for strongly parabolic equations, the result (|9.1U|) is stronger than 
the conclusion of Theorem 12.41 
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The proof of (|9.4j) is based on the analysis of the propagator 



(9.11) 



u a (t) = u I(\a\ = 0) + J (Au a (s) + f(s)I(\a\ = Q)^ds 



+ / Y V a H M kU a -u k) (s) + g k (s)I(\a\ = l))mi(s)ds. 



We consider three particular cases: (1) / = g k = (the homogeneous equation); (2) 
u o = 9k = 0; (3) uq = f = 0. The general case will then follow by linearity and the triangle 
inequality. 

Denote by (&t,s, t > s > 0) the semi-group generated by the operator A(t); := &t,o- One 
of the consequence of Theorem 12.41 is that, under Assumption Al, this semi-group exists 
and is strongly continuous in H. 

Consider the homogeneous equation: / = g k = 0. By Corollary 16.61 

rt «„ r-S2 

(9.12) Yl IM*)llir= / / ■■■/ \\$t,s n M kn ---$ S2 , Sl M kl $ Sl u \\ 2 H ds n , 

\a\=n k u ...,k n >l J ° J ° J ° 

where ds n = ds\ . . . ds n . Define F n (t) = Yl\ a \=n \\ u a(t)\\iii n > 0. Direct application of 
()9.3j) shows that 

(9.13) j t F (t) < C 2 F (t) - \\M k $ t u \\ 2 H . 

k>l 

For n > 1, equality (|9.12|) implies 

"* rSn-l 



(9.14) 
By 



T F n (t)= E / / •••/ \\M kn ^ Sn _ 1 ---M kl ^ Sl u \\ 2 Hds n ' 

rt rs n r-s-2 

+ Y ... (A<S>t,s n Mk n ---$s 1 u 0l <S>t,s n Mk n ---$s 1 u }ds n . 

u. Z m^O Jo Jo 



fcl,...,fe n >l ' 



rt rs n rS2 

V / / .../ {A<S>t, Sn M kn ...$ sl uo,<f>t,s n M kn ...$ sl u }ds n 
k u ...,k n >i Jo Jo Jo 

rt rs n rS2 

(9.15) <- Y J J ■■ ] \\ M k n+ i^t, Sn M kn ...M kl ^ Sl u \\ 2 H ds n 



fci,...,fc n+ i>i 



'0 JO JO 

"t rs n rs2 



rt rs n rs2 

+C 2 Y / / .../ \\^t,s n M kn ...M kl ^ Sl uo\\ 2 H ds n 
,. ^1 Jo Jo Jo 



fel,...,fcn>l 



As a result, for n > 1, 



^nC*) < C 3 F n (t) 



H rs n -\ rS2 

i.. £ ^1 Jo J o Jo 



/•t fSn /-S2 

E //•••/ ll^ fcn+1 $ tlS „/W ifcn .../W fel $ Sl no||l f ^ n . 

. JO JO JO 
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Consequently, 

, N N 

(9.i7) S EE kwiih < E kwiiIt, 

n=0|a|=n n=0\ a \=n 

so that, by the Gronwall inequality, 

N 

(9-18) £ Yl ha(t)\\ 2 H < e C2t \\M 2 H 

n=0 \a\=n 

or 

(9.19) m<t)\?H<^\\M?H' 

Next, let us assume that uq = g k = 0. Then the propagator (|9. llf) becomes 

(9.20) u a (t)= [ {Au a {s) + f(s)l(\a\=0))ds + [ V J^M k u a - m (s)mi(s)ds. 

Jo Jo i,k 

Denote by u^(t) the solution corresponding to a = 0. Note that 

IKo)(*)IIh = 2 [ (Au< )(s),u, 0) (s))ds + 2 I (f(s),u {0) (s))ds 
Jo Jo 

<C 2 f \\ U{0) ( s )f H d S - f^WMkU^WUs + Cf f \\f(s)\\ 2 v ,ds. 
Jo Jo k>1 Jo 

By Corollary 16.61 



(9.21) \\u a (t)f H = E //•••/ \\$t,s n M kn ...M klU{0) ( Sl )f H ds n 

M=n kl ,...,kn>l J ° J ° J ° 

for n > 1. Then, repeating the calculations ()9.14j) -- (|9.16j) . we conclude that 

(9.22) E iK(*)iilr<c/ / ii/(s)ii 2 ^+c 2 / y, E KWiilr*, 

n=l\a\=n J ° J ° n=l\a\=n 

and, by the Gronwal inequality, 

(9.23) n<t)\\H <C f e c * f \\f{s)f v ,ds. 

Jo 

Finally, let us assume that uq = f = 0. Then the propagator (|9.11f) becomes 

u a (t) = J Au a (s)ds 

(9-24) t / \ 

+ J ( E y a i M kU a -(i,k){s) +g k (s)I(\a\ = 1) I rrii{s)ds. 

Even though it a (i) = if a = 0, we have 
(9.25) U(j fc) = / $ t , s gk(s)mi(s)ds, 
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and then the arguments from the proof of Corollary 16.61 apply, resulting in 

"t r s n r s % 



E)H U «(*)I& = E / / •••/ \\®t, Sn M kn ...$ S2 ,s 1 g kl (s 1 )f H ds 

\a\=n k lt ...,k n >l Jo J ° J ° 



for n > 1. Note that 



E WMt)\\ 2 H = E I \\9k{s)\\ 2 H ds + 2Y I {A$ t ,s9k(s),®t,s9k(s))ds. 

\a\=l fc>l^° k>l^° 

Then, repeating the calculations (|9.14j) - (|9.16|) . we conclude that 

(9.26) Y Y KOIIh < E / \\gk(s)f H ds + c 2 / Y Y Ktoll**. 

n=l\a\=n k>l J ° Jo n=l \a\=n 

and, by the Gronwal inequality, 

(9.27) E\\u(t)\\ 2 H <e c ^Y I \\9k(s)f H ds. 

k>i Jo 

To derive (|9~I|) . it remains to combine (|9.19j) . ()9.23|) . and (|9.27|) with the inequality (a + 6 + 
c) 2 <3(a 2 + 6 2 + c 2 ). 

Representation (|9.5j) of the Wiener chaos solution as a sum of iterated ltd integrals now 
follows from Corollary 16.61 Theorem 19. II is proved. □ 

rT 

Corollary 9.2. If / Hiia^Hycfe < do, then sup H^a^lln < °°- 



Proof. The proof of Theorem 19.11 shows that it is enough to consider the homogeneous 
equation. Then by inequalities (|9.15j) - ([9.16|) . 

Y Y im*)iih= E 

(g _ 28) i=n + l\ a \=Z ^=^H 

<e c * T E / / r •■■ n\\M kn+1 $ t ,s n M kn ...$sM\Hds n dt. 

k u ...,k n+1 >l J ° J ° J ° J ° 

By Corollary 16.61 
-T 

' 2 >ds 

*2 



\\u a (s)\\yds 

(9-29) J ° T t Sn 

= Y Y ///•••/ ll^fc^Mn^fcn---^1^0|lv-^ n ^<OO. 

n>lk 1 ,...,k n >l J ° J ° J ° J ° 

As a result, (|93|) and (|9~2U)) imply 

lim / / / ... / \\Mk n+1 ®t,s n Mk n ...M kl ® sl uo\\ 2 H ds n dt = 0, 
n ->°°Jo Jo Jo Jo 

which, by (|9.28|) . implies uniform, with respect to t, convergence of the series ^2 a& j \\u a (t)\\ 2 H . 
Corollary 19.21 is proved. □ 
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Corollary 9.3. Let dij , b{, c, (Tjfc, v k be deterministic measurable functions of (t,x) so that 

\aij(t,x)\ + \bi(t,x)\ + \c(t,x)\ + \a ik (t,x)\ + \v k {t,x)\ < K, 
i,j = l,...,d, k>l, xeM d , 0<t<T; 



i(t,x) - ^a ik (t,x)a jk (t,x)^j yiyj > 0, 



x,y G R d , < t < T; and 

y~] \vk(t,x)\ 2 <C v <oo, 
fe>l 

x G R , < t < T. Consider the equation 

(9.30) du = (Di(aijDju) + b{DiU leal f)dt + (a ik DiU + v k u + g k )dw k . 

Assume that the input data satisfy u G L 2 (R d ), f G L 2 ((0, T); iJ^ 1 (R d )), 
J2 k >i ll5'fclli 2 ((o T)xR d ) < °°' anc ^ fiere exists an e > so that 

a ij (t,x)y i y j > e\y\ 2 , x,y£R d , < t < T. 

Then there exists a unique Wiener Chaos solution u = u(t,x) of 19.,'JU\) . The solution has 
the following regularity: 

(9.31) u(t, •) G L 2 (W; L 2 (R d )), < t < T, 
and 

E l|«lli 2 (Rrf)(*) < C *(lKHi 2(Md) + \\f\\ 2 L2 ^ T) , H -\ 

k>l 

where the positive number C* depends only on C U ,K,T, and e. 
Remark 9.4. 

(1) If holds instead of \9.3\) . then the proof of Theorem \9.1\ in particular, l.9.i.5)) - 

\9.1b]) . shows that the term E||-u(t)||^ in the left-hand- side of inequality \9.4\) can be replaced 
with 

r-t 



(9.32) v - M , 

Il 2 ((o,t)xI 



E( \\u(t)\\ji +£ J \\u(s)\\ 2 v ds 

(2) If f = g k = and the equation is fully degenerate, that is, 2(.A(t)t>, u)+^ fc>1 H-MfcW^II/f = 
; t G [0, T], then it is natural to expect conservation of energy. Once again, analysis of 
\9.15\) - l9.1b}) shows that equality 

E\\u(t)\\ 2 H = \\u \\ 2 H 

holds if and only if 

rT r-t ps n i-S2 

lim / / / ••• / \\M kn+1 $t,s n M kn ...M kl $ Sl uo\\ 2 Hds n dt = 0. 
n ^°°Jo Jo Jo Jo 

The proof of Corollary \9.2\ shows that a sufficient condition for the conservation of energy 
in a fully degenerate homogeneous equation is E \\u(t)\\ydt < oo. 

One of applications of the Wiener Chaos solution is new numerical methods for solving the 
evolution equations. Indeed, an approximation of the solution is obtained by truncating the 
sum Ylaej u a(t)S, a . For the Zakai filtering equation, these numerical methods were studied 
in |25l 1261 12*7] : see also Section fTTI below. The main question in the analysis is the rate of 
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convergence, in n, of the series S n >iS|Q,| =n ||w(t)||#. In general, this convergence can be 
arbitrarily slow. For example, consider the equation 

du = —u xx dt + u x dw(t), t > 0, i£R, 

in the normal triple (H\ (R), L2 0R), H^ 1 ^)), with initial condition u\t=o = u o S Z^QR)- It 
follows from (|9~T2]l that 

j.n /• 

*■»(*) = E IMIW*) = ^ J R \y\ 2n e- y Wdy, 

\a\=n 

where uq is the Fourier transform of u$ . If 

\Mv)\ 2 = {1+ ] yl 2 )r 7>iA 

then the rate of decay of -F n (i) is close to n~( 1+27 ^ 2 . Note that, in this example, 

E IMIi 2 (R)W = INIIi 2{R) - 

An exponential convergence rate that is uniform in H^oIIh ^ s achieved under strong parabol- 
icity condition (|2.5[) . An even faster factorial rate is achieved when the operators A4 k are 
bounded on H. 

Theorem 9.5. Assume that the there exist a positive number e and a real number Co so 
that 

2{A(t)v,v) +Y,\\M k (t)v\\ 2 H + e\\v\\ 2 v <C \\vf H , t € [0,T], v £ V. 

k>l 

Then there exists a positive number b so that, for all t £ [0, T], 



(9-33) \\u a (t)\\ 2 H 



uo\\h 



(1 + 6)" 

|a|=n 

If, in addition, Y^,k>i ll^fe^Mlj? — C3 1 1 V 9 1 1 li" ? then 

(9.34) £ | M 0|&<^e^||«o|&- 

|a|=n 

Proof. If C* is from (|9.6|) and b = e/C* , then the operators Vl + kA/U satisfy 
2(A(t)v, v ) + (l + b)Y^ \\M k (t)f H < C \\vf H . 

k>l 

By Theorem 19.11 

ft fSn fS2 

(l + 6) n V //.../ ||$ Mn M fen ...M fel $ Sl « ||Hds n < INHtf, 

k u ...,k n >l Ja J ° J ° 

and (19331 follows. 

To establish (|93i|l . note that, by (JOJ), 

||**/||!r<e Clt ||/||3!r, 

and therefore the result follows from (|9.12|) . Theorem 19.51 is proved. □ 



The Wiener Chaos solution of ()9.1|) is not, in general, a solution of the equation in the 
sense of Definition 12.21 Indeed, if u L2($l x (0, T); V), then the expressions (Au(s),<p) 
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and (.Mfcu(s), (p)u are not defined. On the other hand, if there is a possibility to move the 
operators A and M. from the solution process u to the test function ip, then equation (|9.1|) 
admits a natural analog of the traditional weak formulation (|2.4|) . 

Theorem 9.6. In addition to A1-A3, assume that there exist operators A*(t), J\A* k {t) and 
a dense subset Vq of the space V so that 

(1) A*(t)(V ) c H, M* k (t)(V ) QH,te [0,T} ; 

(2) for every ve V, if 6 Vq, and t G [0,T], (A(t)v,<p) = (v, A* (t)ip) H , (M k (t)v,<p) H = 
(v,M* k (t)^) H . 

Ifu = u(t) is the Wiener Chaos solution of \9. then, for every ip 6 Vq and every t £ [0, T], 
the equality 

(u(t),tp) H = (u ,ip) H + / (u(s),A*(s)ip) H ds + / (f(s),cp)ds 
(9.35) Jo Jo 

+ / (u(s),^fc(s)v3) J ffdw fc (s) + / {gk(s),<p)Hdw k (s) 



o JO 



ZioWs in L 2 (W). 

Proof. The arguments are identical to the proof of Theorem I8.6f 2). □ 



As was mentioned earlier, the Wiener Chaos solution can be constructed for anticipating 
equations, that is, equations with J^-measurable input data. With obvious modifications, 
inequality (|9.4|) holds if each of the input functions uq, f, and g k in (J9.1|) is a finite linear 
combination of the basis elements £ a . The following example demonstrates that inequality 
(J9.4j) is impossible for general anticipating equation. 

Example 9.7. Let u = u(t,x) be a Wiener Chaos solution of an ordinary differential 
equation 

(9.36) du = udw(t),0 < t < 1, 

with uq = Yla&j a a£,a- For n > 0, denote by (n) the multi-index with ot\ = n and aj = 0, 
i > 2, and assume that a( n ) > 0, n > 0. Then 

(9.37) E« 2 (l)>C^e^. 

n>0 

Indeed, the first column of propagator for a = (n) is Um\(t) = a( ) and 

u (n) (*) = O(n) +Vn Ut n -l)(s)ds, 

Jo 

so that 



(n) 



, Q y/(n-k)\k\ 
Then^ n) (l)>ELo ©^and 



n>0 n>0 \fc>0 v 



1 fn + k\ \ 2 

a (n)- 
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Since 



the result follows. 



k>0 



k\\ n 



k>o v 



The consequence of Example 19.71 is that it is possible, in (|9.1j) . to have uq £ L^W;-?/) for 
every n, and still get E||n(t)||^ = +00 for all t > 0. More generally, the solution operator 
for (|9,1|) is not bounded on any Z/2,Q or (5)_ Pi _ 7 . On the other hand, the following result 
holds. 

Theorem 9.8. In addition to Assumptions Al, A2, let u$ be an element of T>'(W;H), 
f, an element of V'(W; L 2 ((0, T), V' )), and each g k , an element of V'(W; L 2 ((0, T), H)). 
Then the Wiener Chaos solution of equation i9.1\) satisfies 



(9.38) 



£ \\Ua(t)\\H 



\U0a\\H + 



t \ 1/2 

ll/«(*)llv"<k 



where C > depends only on T and the numbers 5, C\, and C2 from and \9. 

Proof. To simplify the presentation, assume that / = g k = 0. For fixed 7 G J, denote by 
u(t; cp; 7) the Wiener Chaos solution of the equation (|9.1|) with initial condition u(0; (p; 7) = 
Denote by (0) the zero multi-index. The structure of the propagator implies the 
following relation: 



(9.39) 



u, 



K+7 (i;^ 7 ) _ u " ( t; v% ; ^) 



V(a + 7)! 
Clearly, u Q (i; 92; 7) = if |a| < | -y | . If 



'al 



(S)-i,o(H) 



a6j 



" 1 ' 



al 



then, by linearity and triangle inequality, 



|u(t)[| 



(5)-i,o(H) 



< 



II«(*;«07;7)II(s)_i,o(h)- 



7SJ 



We also have by (|9.39f) and Theorem 19.11 

lk(*;«0 7 ;7)ll(5)_ 1 , (fl-) 

< E 



«U;^X;(0) 



« 1*5^5(0) 



'7! 



(S)-i,o(H) 
2 



< e 



C 2 i 



II«0tIIh 



7! 



Inequality (j9,38|) then follows. Theorem 19.81 is proved. 



□ 



40 S. V. LOTOTSKY AND B. L. ROZOVSKII 

Remark 9.9. Using Proposition \ 7. 6 <\ and the Cauchy- Schwartz inequality, iy.3H\) can be 

re-written in a slightly weaker form to reveal continuity of the solution operator for equation 
19. 1\) from (<S)_i )7 to (5)_i f o for every 7 > 1: 

ll u (*)ll?S)-i,o(fl) ^ C [\\M\\s)-^{H) + J o ll/( s )ll(5)-i, 7 (V') ds 

+ Y1 I \\9h{s)\\ 2 {s) {H) ds\. 
k>l Jo J 

10. Probabilistic Representation of Wiener Chaos Solutions 

The general discussion so far has been dealing with the abstract evolution equation 

du = (Au + f)dt + y~^(.Mfcn + gk)dwk- 

k>l 

By further specifying the operators A and Mk, as we ll as the input data uq, f, and gk, it is 
possible to get additional information about the Wiener Chaos solution of the equation. 

Definition 10.1. For r E R, the space ^2,(r) = -^2,(r)0^ d ) ^ s ^ e collection of real-valued 
measurable functions so that f £ L 2 m if and only if J Rd |/(x)| 2 (l + |x| 2 ) r <ix < 00. The space 
^ = ^ (M. d ) is the collection of real-valued measurable functions so that f £ , % if 
and only if f and all the first-order generalized derivatives Dif of f belong to £2,(7-) ■ 

It is known, for example, from Theorem 3.4.7 in 02j, that -^2,(r) is a Hilbert space with 
norm 

ll/llo,(r)= / \f(x)\ 2 (l + \x\ 2 Tdx, 
and H}, is a Hilbert space with norm 

d 

ll/lll,(r) = ll/llo,(r)+X;ilA/||o,(r)- 
i=l 

Denote by ^ ne dual of ^2 (r) w *th respect to the inner product in L 2 ^ r y Then 

(i?2 r r y ^2,(r)i ^2(r0 * s a norma l triple of Hilbert spaces. 

Let F = (f2, J 7 , {Ji}{>o, P) be a stochastic basis with the usual assumptions and = 
Wk{t), k > 1, i > 0, a collection of standard Wiener processes on F. Consider the linear 
equation 

(10.1) du = (aijDiDjU + biDiU + cu + / )<it + {a^BiU + z^u + gk)dwk 

under the following assumptions: 

BO All coefficients, free terms, and the initial condition are non-random. 

Bl The functions ay = a,ij(t,x) and their first-order derivatives with respect to x are 
uniformly bounded in (t, x), and the matrix (aij) is uniformly positive definite, that 
is, there exists a 5 > so that, for all vectors y £ M. d and all (t, 2), a^yiyj > 5\y\ 2 . 

B2 The functions = bi(t,x), c = c(t,x), and I/& = Uk(t,x) are measurable and 
bounded in (t, x). 

B3 The functions a^. = o~ik(t,x) are continuous and bounded in (t, x). 
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B4 The functions / = f(t,x) and gk = gk(t,x) belong to L2((0,T); L 2 ,(r)) for some 
rEl. 

B5 The initial condition uq = uq{x) belongs to L^At)- 

Under Assumptions B2-B4, there exists a sequence Q = {q k , k > 1} of positive numbers 
with the following properties: 

PI The matrix A with Ay = a y - — (1/2) ^ fc>1 qk&ikOjk satisfies 

Aj(t,x)yiVj > 0, 

x,y G M d , < t < T. 
P2 There exists a number C > so that 

^ fsup|g fc ^ fc (t,x)| 2 + / \\qk9k\\ P ^ r) (t)dt) < C. 
fc > 1 V t,x Jo ) 

For the matrix A and each i, x, we have Ay(t, x) = aik(t,x)ajk(t,x), where the functions 
&ik are bounded. This representation might not be unique; see, for example, Theorem 
III. 2. 2] or |441 Lemma 5.2.1]. Given any such representation of A, consider the following 
backward Ito equation 



(10.2) 



X t , x ,i (s) = Xi+ Bi (r, X t , x (r)) + Qk^ik (r, X t>x (r)) dw k (r) 

•* s k>l 

+ I a ik (r, X tiX (r)) cM> fc (r) ; s G (0, t), t G (0, T], t - fixed, 

JS 



where Bi = bi — X^fc>i 9fc cr ifc I/ fc and & > 1, are independent standard Wiener processes 
on F that are independent of Wk, k > 1. This equation might not have a strong solution, 
but does have weak, or martingale, solutions due to Assumptions B1-B3 and properties 
PI and P2 of the sequence Q; this weak solution is unique in the sense of probability law 
HH Theorem 7.2.1]. 

The following result is a variation of Theorem 4. 1 in |29j . 

Theorem 10.2. Under assumptions BO B5 equation \10. 1\) has a unique 

w(i?2 ( r )i ^2~(r)) Wiener Chaos solution. If Q is a sequence with properties PI and P2, 

then the solution of belongs to 



L 2 ,Q (W; L 2 ((0, T); H^)) f| L 2 , Q (W; C((0, T); L 2 , (r) )) 
and has the following representation: 



u(t,x) = Q- 1 E[ I f(s,X t<x (s)h(t,s,x)ds 

(10.3) 



k>l 



+ / 9fc5fe(s, A tj:r (s))7(t,s,x)du' fc (s) + uo(Ai ia; (0))7(t,0,a;) 



, i<T, 
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where X t:X {s) is a weak solution of \1(). and 



t i 

q k u k (r, X t ^ x (T))dw k (T) 



j(t,s,x) = exp [ / c(T,X tjX (T))dT + ^2 f 

-lj^ql\Mr,XUr))\ 2 dr\ 



Proof. It is enough to establish (|1U.3|) when t = T. Consider the equation 

(10.5) dU = (aijDiDjU + b^U + cU + f)dt + ^{a ik D t U + v k U + g k )q k dw k 



k>l 



with initial condition U(0, x) = Uq(x). Applying Theorem 12.41 in the normal triple 
(H^ ( r yL2 i ( r ),H~}s), we conclude that there is a unique solution 

U e L 2 (W;L 2 ((0,T); J ff 2 1 (r) )) f|L 2 (W; C((0, T); L 2 , (r) )) 

of this equation. By Proposition 17.41 the process u = Q~ 1 U is the corresponding Wiener 
Chaos solution of (|10.1|) . To establish representation ()10.3() . consider the S-transform Uh 
of U. According to Theorem 18. 11 the function Uh is the unique w(H^ ,y solution of 

the equation 

(10.6) dU h = (aijDiDjUh + biDiU h + cU h + f)dt + ^2(a ik DiU h + v k U h + gk)qkh k dt 

k>i 

with initial condition Uh\t=o = uo- We also define 



(10.7) 



Y(T,x)= f T f(s,X T . x (s)h(T,s,x)ds 
Jo 

+ 9k(s,X T!X (s))j(T,s)q k dw k (s)+u (X TtX (0))^(T,Q,x). 



fc>i ■ 

By direct computation, 

E(E (S(h)Y(T,x)\J^)) = E (£(h)Y(T, x)) = E'Y(T,x), 

where E' is the expectation with respect to the measure d¥' T = £(h)d¥T and Pt is the 
restriction of P to T^f . 

To proceed, let us first assume that the input data uq, f, and g k are all smooth functions with 
compact support. Then, applying the Feynmann-Kac formula to the solution of equation 
(|10,6[) and using the Girsanov theorem (see, for example, Theorems 3.5.1 and 5.7.6 in |15j). 
we conclude that Uh(T,x) = K'Y(T,x) or 

E (£(h)EY(t, x)\J^) = E {£ (h) U(T, x)) . 

By Remark 17.81 the last equality implies U(T,-) = E (Y(T, Ol-^r ) as elements of 
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To remove the additional smoothness assumption on the input data, let Uq, f n , and g% be 
sequences of smooth compactly supported functions so that 



(10. 



hm mk - u%\\i 2Mm + y o ||/ - f n \\i 2Xr)m (t)dt 
+ E / vl\\9k-gnl Mr)( R d) (t)dt) =o. 

k>l J0 I 



Denote by U n and Y n the corresponding objects defined by Q10.5JI and (|10.7|) respectively. 
By Theorem l9.il we have 



(10.9) 



lim E\\U - U n \\ 2 r 



(T) = 0. 



To complete the proof, it remains to show that 



(10.10) 



lim E 

n^oo 



E (Y(T,-) - Y n (T, •) 



-r-W 
J- T 



0. 



-2,(r) 



To this end, introduce a new probability measure P T by 

cT 



dF? = exp 



v k (s,X% x (s))q k dw k (s) 



2 r T,^(^ x ?A s ^ 2ds ) dF T. 
j ° k >i j 



By Girsanov's theorem, equation (|10.2|) can be rewritten as 



X T , x ,i (s) =Xi+ ^2 a ik (r, X T , X (t)) h k (t) q k dr 
Js k>l 

(10.11) + f (bi + ^2q 2 k (Ti k u k )(T,X T ^(T))dT 

Js k>i 

+ f ^2qkcrik(r,X TtX (T))dw'^.(T) + f d ik (t, X t , x (t)) dw" k (r) 



k>l 



where w'l and w" k are independent Winer processes with respect to the measure P^. Denote 
by p(s,y\x) the corresponding distribution density of Xt, x (s) and write £{x) = (1 + |x| 2 ) r . 
It then follows by the Holder and Jensen inequalities that 

2 



(10.12) 



E 



E 



(j\\T,s,-)(f - r)(s,X T ,(s))ds 



J- T 







E ( 7 2 (r,s,x)(/ - p) 2 ( S ,X T , x (s))) ds £(x)dx 



< K 2 

■T 



- f n )\s,X T>x {s))ds) i{x)dx 







(f(s,y) - f n (s,y)) 2 p(s,y\x)dy ds £(x)dx, 



where the number K\ depends only on T, and the number K<i depends only on T and 
sup( ti2 .) |c(i,x)| + Ylk>i Ik su P(t,x) \ v k{t,x)\ 2 . Assumptions B0-B2 imply that there exist 
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positive numbers K3 and K4 so that 

(10.13) P(s,y\x)< 



K 3 I \x-y\ 

exp — A4- 



(T - s) d / 2 



T-s 



see, for example, jB]. As a result, 



and 



r 



JM. d 



p(s,y\x)£(x)dx < K 5 £(y), 



(f(s, y) ~ f n (s, y)) 2 p(s, y\x)dy ds £(x)dx 



<k 5 J q ||/ -rill. 



(10-14) 

11 f f n \ 

where the number K$ depends only on K^,K^, T, and r 
Calculations similar to (|10.12|) - (|l().14jl show that 

2 

L 

(■T 



(s)ds — > 0, n — > 00, 



E 



E 7 2 (r,0,0(no-0(X T) .(0)) 



(10.15) 



+ E 



E 



^2(9k ~ 9k)(s, X T: .(s)h{t, s, -)q k dw k (s) 



k>l 







as n — > 00. Then convergence ()1().1()|) follows, which, together with (|l().9j) . implies that 
f7 (T, •) = E (C/ Q (T, -)|J^f ) as elements of L 2 (O; L 2j(r) (IR d )) . It remains to note that u = 
Q~ l U . Theorem 110.21 is proved. □ 



Given / £ ^2,(r)> we sa Y t na t / > if and only if 

/ f(x)cp(x)dx > 

for every non-negative if 6 C§°(]R rf ). Then Theorem 110.21 implies the following result. 

Corollary 10.3. In addition to Assumptions BO B5, let uq > 0, / > 0, and g k = /or 
a// k > 1. T/ien u > 0. 



□ 



Proof. This follows from (|10.H|) and Proposition 17. 1 ll 

Example 10.4. (Krylov-Veretennikov formula) 
Consider the equation 

d 

(10.16) du = (aijDiDjU + biDiu) dt + aikDiudw k , u (0, :r) = no (x) . 

k=l 

Assume B0-B5 and suppose that aij(t,x) = ^aik(t,x)ajk(t, x). By Theorem 19.11 equation 
(|10.16j) has a unique Wiener chaos solution so that 



Ellnl 



L-2 



(t)<C*||« |li 
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and 



x d "t r s n fS 2 



(10.17) 



■(t,x)= ^2 e u a (t,x)€ a = u (x) +e ^2 / / ■ • •/ 

n=l\ a \=n n=l k lt ...,k n =l J ° J ° " 

®t,s n Vjk n Dj ■ ■ ■ $ S2>sl a ikl D i § slt0 uo(x)dw kl (s 1 ) ■ ■ ■ dw kn (s 



where 3>t )S is the semi-group generated by the operator A = aijDiDjii + biD^u. On the 
other hand, in this case, Theorem 1 1 U . 21 yields 

u(t,x)=E(u (X t>x (0))\^y 

where W = (w%, ...,w^) and 

rt d < 

X t , x ,i (s) =Xi+ hi (r, X t>x (r)) dr + ^ a ik (r, X t>x (r)) dw k (r) 

J S ,.. i 



(10.18) 



fc=i 

s G (0,t), i G (0,T], i- fixed. 



Thus, we have arrived at the Krylov-Veretennikov formula |201 Theorem 4] 

°° d „t „ Sn „ S2 

E( uo (x t , x (o))\ry) = uo(x) + J2 £ //•••/ 

( 10 - 19 ) n=lk u ...,kn=l Jo J ° J ° 

^t,s n crj kn Dj • • • <& S2jSl (Ti kl Di$ Slfi u (x)dw kl (si) • • • dw kn (s n ). 



11. Wiener Chaos and Nonlinear Filtering 

In this section, we discuss some applications of the Wiener Chaos expansion to numeri- 
cal solution of the nonlinear filtering problem for diffusion processes; the presentation is 
essentially based on 

Let (p., T , P) be a complete probability space with independent standard Wiener processes 
W = W(t) and V = V(t) of dimensions d\ and r respectively. Let Xq be a random variable 
independent of W and V. In the diffusion filtering model, the unobserved d - dimensional 
state (or signal) process X = X(t) and the r-dimensional observation process Y = Y(t) are 
defined by the stochastic ordinary differential equations 

dX(t) = b(X(t))dt + a(X(t))dW(t) + p{X(t))dV(t), 
(11.1) dY(t) = h(X(t))dt + dV(t), < t < T; 

X(0) = X , Y(0) = 0, 

where b(x) G R d , a(x) G R dx *, p(x) G R dxr , ft(a?) G W. 



Denote by C n (R ) the Banach space of bounded, n times continuously differentiable func- 
tions on R rf with finite norm 

C n (Rd) = sup \f(x)\ + max sup \D k f(x)\. 

xm d l<k<n xi 



Assumption Rl. The the components of the functions a and p are in C 2 (R d ), the com- 
ponents of the functions b are in C 1 (R), the components of the function h are bounded 
measurable, and the random variable Xq has a density uq. 
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Assumption R2. The matrix aa* is uniformly positive definite: there exists an e > so 
that 

d d\ 

^2 ^2vik(x)vjk{x)yiyj > e\y\ 2 , x,y G 

i,j=l k=l 



t,d 



Under Assumption Rl system has a unique strong solution |15[ Theorems 5.2.5 and 

5.2.9]. Extra smoothness of the coefficients in assumption Rl insure the existence of a 
convenient representation of the optimal filter. 

If / = fi x ) is a scalar measurable function on W d so that 

supQ <t<T K\f(X(t))\ 2 < co, then the filtering problem for Ijll.lj) is to find the best mean 
square estimate ft of f(X(t)), t <T, given the observations Y(s), < s <t. 

Denote by Tj the <7-algebra generated by Y(s), < s < t. Then the properties of the 
conditional expectation imply that the solution of the filtering problem is 

f t = E(f(X(t))\Ff). 

To derive an alternative representation of ft, some additional constructions will be necessary. 
Define a new probability measure P on (O, J-) as follows: for A G T , 

P(A) = / Z T l a¥, 
J A 

where 

rt i rt 



exp /I* h*{X{s))dY{s) -I J* \h(X(s))\ 2 ds\ 



(here and below, if £ G M fc , then £ is a column vector, = (Ci, . . . ,Cfc); an d |C| 2 = C*0- 
If the function h is bounded, then the measures P and P are equivalent. The expectation 
with respect to the measure P will be denoted by E. 

The following properties of the measure P are well known |14l I42j : 

PI. Under the measure P, the distributions of the Wiener process W and the random 
variable Xq are unchanged, the observation process Y is a standard Wiener process, 
and, for < t < T, the state process X satisfies 

dX{t) = b(X{t))dt + a(X(t))dW(t) + p(X(t)) {dY{t) - h{X(t))dt) , 
X(0) = X ; 

P2. Under the measure P, the Wiener processes W and Y and the random variable Xq 

are independent of one another; 
P3. The optimal filter ft satisfies 



(11-2) /, 



E [f(X(t))Z t \f t 



E[Z t \F} 



Because of property P2 of the measure P the filtering problem will be studied on the 
probability space (£l,J-,F). In particular, we will consider the stochastic basis F = 
{Q,T,{T Y } <t<T,F} and the Wiener Chaos space L2(Y) of ^"^-measurable random vari- 
ables rj with E|?7| 2 < oo. 
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If the function h is bounded, then, by the Cauchy-Schwarz inequality, 



(11.3) E\rj\ < C(h,T)y/E\ri\ 2 , rj G L 2 (Y). 

Next, consider the partial differential operators 

Cg{x)= l -Y, ((^K(^ + W^(*))«)^ + EM»)^i 

i,j=l % ^ 8=1 1 



Mig(x) = hi(x)g(x) + 2J ft; 0*0 g > / = 1, . . . , r; 

i=l 

and their adjoints 



1 d d 2 

£*9{x) = nYl a a (M^V*^))^^) + (j>(?)p*(x))ijg(x)) 



dxidxn 

i,j=l ■> 





(M^M^)); 



i=i 9 ^ 



M*[g(x) = hi(x)g(x) - V — (p u (x)g(x)) , / = 1, . . . , r. 

i=i 

Note that, under the assumptions Rl and R2, the operators £, £* are bounded from H\ 
to H 2 1 (R d ), operators M,M* are bounded from H%(R d ) to L 2 (IR d ), and 

r 

(n.4) 2<r«,«) + ^11^11^) + Eiiuii^^j <c|Ml| 2(Rd) , ^^rt, 

i=i 

where (•, •) is the duality between H^Mf 1 and -ff^" 1 (IR d ). The following result is well known 
|H Theorem 6.2.1]. 

Proposition 11.1. In addition to Assumptions Rl and Rl suppose that the initial density 
Uq belongs to L^IR 01 ). Then there exists a random field u = u(t,x), t G [0, T], x G M. d , with 
the following properties: 

l.ue L 2 (Y; L 2 ((0, T); H$(R d ))) n L 2 (Y; C([0, T], L 2 (M d ))). 

12. T/ie function u(t, x) is a traditional solution of the stochastic partial differential equation 



du(t,x) = £*u(t,x)dt + ^^M*u(t,x)dYi(t), 0<t<T, xG 
^ " ' z=i 
w(0, x) = no(x). 

3. The equality 

(11.6) I [f(X(t))Z t \J* ] = f f(x)u(t,x)dx 

holds for all bounded measurable functions f . 
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The random field u = u(t, x) is called the unnormalized filtering density (UFD) and the 
random variable <fit[f] = E [f(X(t))Zt\J-Y] , the unnormalized optimal filter. 

A number of authors studied the nonlinear filtering problem using the multiple ltd integral 
version of the Wiener chaos [21l21 [ l3*n il46| etc.]. In what follows, we construct approximations 
of u and 4>t[f] using the Cameron-Martin version. 

By Theorem 

(11.7) u(t,x) = ^2 u a (t,x)£ a , 
where 

(11.8) e a = -^=TT £«= / mi (t)dY k (t), k = l,...,r; 

as before, H n (-) is the Hermite polynomial (|3.3|) and m; S m, an orthonormal basis in 
L2((0,T)). The functions u a satisfy the corresponding propagator 

d 

—u a (t,x) = C*u a (t,x) 

(11.9) +^^M* k u a - m (t,x)mi(t), < t < T, x£R d ; 

k.i 

u(0,x) = uo(x)I(\a\ = 0). 

Writing 

fa(t) = / f(x)u Q (t,x)dx, 
JR d 

we also get a Wiener chaos expansion for the unnormalized optimal filter: 
(11-10) = t€[0,T\. 

For a positive integer N, define 

(11.11) u N (t,x)= ^2 u a (t,x)£, a . 

\a\<N 

Theorem 11.2. Under Assumptions Rl and R2, there exists a positive number v, depend- 
ing only on the functions h and p, so that 

[| Uo ||2 

(H-12) E\\u-u N \\l 2m (t) < ^J^ff , t G [0, T]. 



7/j in addition, p = 0, i/ien i/iere exists a real number C, depending only on the functions b 
and a, so that 



(Ah t) N+1 

(11.13) E^-Uivll^ti) < (JV T +1) , e^NII^,,- / e to.:/"]. 

where hoc = maxfc=i,...,r sup x |/ifc(x) 



For positive integers iV, n, define a set of multi-indices 

J# = {a = (af , fc = 1, . . . , r, i = 1, . . . , n) : \a\ < N}. 
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and let 

(11.14) u%{t,x) = Mt,x)Za- 

Unlike Theorem 111.21 to compute the approximation error in this case we need to choose 
a special basis m — to do the error analysis for the Fourier approximation in time. We 
also need extra regularity of the coefficients in the state and observation equations — to 
have the semi-group generated by the operator C* continuous not only in L2(M d ) but also 
m H%(R d ). The resulting error bound is presented below; the proof can be found in |25j . 

Theorem 11.3. Assume that 

(1) The basis m is the Fourier cosine basis 

(11.15) m 1 (s) = ^=; m*(*) = y| cos ~ , k > 1; < t < T, 



(2) The components of the functions a are in C 4 (M. d ), the components of the functions b 
are in C 3 (M), the components of the function h are in C 2 (R d ); p = 0; uq 6 iTf (R ). 

Then there exist a positive number B\ and a real number Bi, both depending only on the 
functions b and a so that 

~ / (4h T) N+1 T 3 

(11.16) E\\u - u n N \\l 2m (T)<B ie B > T + e W u °Wl 2 m + — Nll^ 

where /i M = max^ =li r sup x \ hk(x)\. 

12. Passive Scalar in a Gaussian Field 

This section presents the results from |29| and [2E] about the stochastic transport equation. 

The following viscous transport equation is used to describe time evolution of a scalar 
quantity 9 in a given velocity field v: 

(12.1) 6(t, x) = vA6(t, x) - v(t, x) ■ V0(t, x) + f(t, x); x£R d , d>l. 
The scalar 9 is called passive because it does not affect the velocity field v. 

We assume that v = v(t, x) G W 1 is an isotropic Gaussian vector field with zero mean and 
covariance 

E(v% x)v\s,y)) = 8(t - s)CV(x - y), 
where C = {C 13 (x),i,j = 1, ...,d) is a matrix-valued function so that C(0) is a scalar 
matrix; with no loss of generality we will assume that C(0) = I, the identity matrix. 

It is known from |221 Section 10.1] that, for an isotropic Gaussian vector field, the Fourier 
transform C = C(z) of the function C = C{x) is 

/ 100 n Af \ A o ( W* . b ( yy T 

(12.2) C(y) = - | a—, + — — / 



(1 + |y|2)(d+a)/2 ^ | y |2 ' d - 1 \ \y\ 2 

where y* is the row vector (yi, . . . ,yd), y is the corresponding column vector, \y\ 2 = y*y; 
7 > 0, a > 0, b > 0, Ao > are real numbers. Similar to |22j . we assume that < 7 < 2. 
This range of values of 7 corresponds to a turbulent velocity field v, also known as the 
generalized Kraichnan model 0; the original Kraichnan model |18j corresponds to a = 0. 
For small x, the asymptotics of C^(x) is (Sij — c lJ |x| 7 ) [HI Section 10.2]. 
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By direct computation (cf. P), the vector field v = (v 1 , . . . , v d ) can be written as 

(12.3) v%x) = ai(x)w k {t), 

where {a k , k > 1} is an orthonormal basis in the space He, the reproducing kernel Hilbert 
space corresponding to the kernel function C. It is known from [22] that He is all or part 
of the Sobolev space H^ d+ ^/ 2 (R d ; R d ). 

If a > and b > 0, then the matrix C is invertible and 

H C = {/ € K d : ^ r{y)C-\y)f{y)dy < ooj = ff( d +7)/ 2 (R d ; M rf ), 

because ||C(y)|| ~ (1 + |y| 2 )-^/ 2 . 
If a > and 6 = 0, then 

ffc = {/et d : J Rd |/(y)| 2 (l + |y| 2 ) (d+7)/2 ^ < oo ; yy *f{ y ) = \y\ 2 f(y)] , 

the subset of gradient fields in H ( d +T)/ 2 (R d ; M rf ), that is, vector fields / for which f(y) = 
yF(y) for some scalar F G #( d +T+ 2 )/ 2 (IR d ). 

If a = and 6 > 0, then 

ff C = {/£K d : £ |/(tf)| 2 (l + \y\ 2 )^ )/2 dy < oo; y*/(l/) = o} , 
the subset of divergence- free fields in H^ d+ ^/ 2 (R d :R d ). 

By the embedding theorems, each erf is a bounded continuous function on M. d ; in fact, every 
a\ is Holder continuous of order 7/2. In addition, being an element of the corresponding 
space He, each a k is a gradient field if b = and is divergence free if a = 0. 

Equation 1)12,1(1 becomes 

(12.4) d6{t, x) = (vA6(t, x) + f{t, x))dt - ^ a k (x) ■ V6(t, x)dw k {t). 

k 

We summarize the above constructions in the following assumptions: 

51 There is a fixed stochastic basis F = (O, T, {^}t>o, P) with the usual assumptions 
and (wk(t),k > 1, t > 0) is a collection of independent standard Wiener processes 
on F. 

52 For each /c, the vector field a k is an element of the Sobolev space 
H ( 2 d+l)/2 (R d ;R d ), < 7 < 2, d > 2. 

53 For all x,y in M d , J2k a k( x ) a l(y) = ^(x — y) so that the matrix- valued function 
C = C(x) satisfies (TTOl and C(0) = J. 

54 The input data 60, f are deterministic and satisfy 

60 G L 2 (K d ), / G L 2 ((0,r);Hj 1 (R <, )) ; 
1/ > is a real number. 
Theorem 12.1. Let Q be a sequence with q k = q < y2u, k > 1. 
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Under assumptions S1-S4, there exits a unique w(H2(M d ), H^ 1 (M d )) Wiener Chaos solu- 
tion of \12.J$ . This solution is an -adapted process and satisfies 

\\a\\2 , \\a\\2 

H t 'llL 2i Q(W;ia((0,T);fri(]R«i))) t" ll p ll£ 2]Q (W;C((0,T);L 2 (]R<*))) 



< C(u,q,T) (||0o|li 2(R d) + ll/lli 2 (( ,T);H 2 - 1 (R<i))) • 



Theorem 112.11 provides new information about the solution of equation 1)12.1)) for all values 
of v > 0. Indeed, if y2u > 1, then q > 1 is an admissible choice of the weights, and, by 
Proposition I7.4f 1). the solution 9 has Malliavin derivatives of every order. If y2u < 1, then 
equation (|12.4|) does not have a square-integrable solution. 

Note that if the weight is chosen so that q = \f2v, then equation (|12.1f) can still be analyzed 
using Theorem EH1 in the normal triple (H^(R d ), L 2 (R d ), H^ 1 ^)). 

If v = 0, equation 1)12.4)1 must be interpreted in the sense of Stratonovich: 

(12.5) du(t, x) = f(t, x)dt - a k {x) ■ V0(i, x) o dw k (t). 

To simplify the presentation, we assume that / = 0. If 1)12.2)1 holds with a = 0, then each 
<jfc is divergence free and 1)12. 5|) has an equivalent ltd form 

(12.6) d8(t, x) = -A9(t, x)dt - a{{x)D % 9{t, x)dw k (t). 

Equation 1)12.6)1 is a model of non- viscous turbulent transport [5]. The propagator for 1)12.6)1 

is 

(12.7) ^e a (t,x) = ±A0 a (t,x) - £ \f^AD j 9 a - {i!k) {t,x)m i {t), 0<t<T, 

i,k 

with initial condition 9 a (0,x) = do(x)I(\a\ = 0). 

The following result about solvability of 1)12.6)1 is proved in 29\ and, in a slightly weaker 
form, in 2i 



Theorem 12.2. In addition to S1-S4, assume that each a k is divergence free. Then there 
exits a unique w(H%(R d ), H^ 1 ^)) Wiener Chaos solution 6 = 9(t,x) of fTOj) . This 
solution has the following properties: 



(A) For every <p G Cg°(M d ) and all t £ [Q,T], the equality 
(12.8) (9, <p)(t) = (9 , <p) + l ! (9, A<p)(s)ds + f\d, alD iV )dw k {s) 



2„ 

holds in h%{T^'), where (•, •) is the inner product in Li2(M. d ). 

(B) If X = Xt >x is a weak solution of 

(12.9) X t>x = x + [ a k (X s>x ) dw k (s) , 

JO 

then, for each t G [0, T], 

(12.10) 9(t,x)=E(9 (X tjX )\^). 

( C) For 1 < p < 00 and r£l, define L p / r \(M. d ) as the Banach space of measurable functions 
with norm 

w^ = I \mni+\x\*Yi*dx 
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is finite. Then there exits a number K depending only on p, r so that, for each t > 0, 

( 12 -n) n0\\l M (t)<e Kt \\e o E 



In particular, if r = 0, then K = 0. 

It follows that, for all s, t and almost all x, y, 

E9(t,x) = 9 a (t, x) I\ a \ =0 
and 

E9(t,x)9(s,y) = ^9 a 9 « (^y)- 

If the initial condition 6 belongs to L 2 (R d ) n L p (R d ) for p > 3, then, by (|12.11|) . higher 
order moments of 9 exist. To obtain the expressions of the higher-order moments in terms 
of the coefficients 9 a , we need some auxiliary constructions. 

For a, € define a + /3 as the multi-index with components Similarly, we define 

the multi-indices \a — j3\ and a A = min(a, 0). We write < a if and only if < for 
all i, A; > 1. If /3 < a, we define 

cA rr a?! 



^ /??!(<*?-/?)!' 

Definition 12.3. We say that a triple of multi-indices (a, /?, 7) is complete and write 
(a, /3, 7) £ A i/ all £/ie entries of the multi-index a + + 7 are et>en numbers and 
\& — 0\ < 7 < a + /3. -For /ixed a, £ J, we write 

A (a) :={7,fi6j: (a, 7, //) G A} 

and 

A(a,/3) :={ 7 G l 7:(a,/3,7)GA}. 

For (a,/3, 7) G A, we define 

'a — + 7\ I //3-a + 7\ 1 / a + /? - 7 



(12.12) ¥(a,/3, 7 ) :=v^w(( a ^ + 7 ) ! (- 



Note that the triple (a, 0, 7) is complete if and only if any permutation of the triple (a, 0, 7) 
is complete. Similarly, the value of \& (a,/3,7) is invariant under permutation of the argu- 
ments. 



We also define 

'7 + ~ 2/A h\ (0 
7 - V ) W 

It is readily checked that if / is a function on J ', then for j,0 & J, 



(12.13) C( 7 ,0,fi) :~- 



1/2 

, /x< 7 A/3. 



(12.14) Yl C(7,Ap)/(7 + /?-2M)= J] /(M)*(7,AA*) 

The next theorem presents the formulas for the third and fourth moments of the solution 
of equation (|12.6|) in terms of the coefficients 8 a . 
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Theorem 12.4. In addition to SI S4, assume that each is divergence free and the 
initial condition Oq belongs to L2(M. d ) DL^W*). Then 

(12.15) E9(t,x)9(t',x')9(s,y) = ^ (a, /?, 7 ) 9 a (t, x)9 (t', x')9^ (s, y) 

(a,/3, 7 )6A 

and 

(12.16) E0(t, x)9(t', x')9 (s, y) 9 (s 1 , y') 

y(a,p,p)y(p, 1 , K )9 a (t,x)9p(t',x')9 1 (s,y)9 K {s',y'). 

p6A(a,/3)nA(7,K) 

Proof. It is known 30j that 

(12.17) = J2 0^,(3,^)^-2^. 



*(«,/?, 7), («,/3,7)eA; 

otherwise. 



Let us consider the triple product £ a £/3£ 7 . By (|12,17|) . 

(12.18) Ed a ^ = E C7^*(«»/3.^) = |* 

Equality (|12.15|) now follows. 

To compute the fourth moment, note that 

p<aA/3 

(12.19) ~ 

= 2J C ( a ^^) z2 C(« + /5-2^,7,P)^+/3+7-2/ 1 -2p- 

A*<aA/3 p<(a+/3-2p)A7 

Repeated applications of (|12,14|) yield 

£«Cj86y = Yl °( a ^^) J2 ^{a + f3-2p, 1 ,p) 

p<aA/3 peA(a+/3-2p,7) 
pSA(a,/3) peA(p, 7 ) 

Thus, 

MSA(a >( 3) pSAO, 7 ) 

J] *(a,/3,p)*(/9,7,K). 

p£A(a,/3)nA(7,K) 

Equality (j!2.1fi[) now follows. □ 



In the same way, one can get formulas for fifth- and higher-order moments. 

Remark 12.5. Expressions h!2.15\) and \12.1b\) do not depend on the structure of equation 
\12.b}) and can be used to compute the third and fourth moments of any random field with 
a known Cameron- Martin expansion. The interested reader should keep in mind that the 
formulas for the moments of orders higher then two should be interpreted with care. In fact, 
they represent the pseudo-moments (for detail see |35j ). 
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We now return to the analysis of the passive scalar equation Q12.4|) . By reducing the 
smoothness assumptions on o~ k , it is possible to consider velocity fields v that are more 
turbulent than in the Kraichnan model, for example, 

(12.20) v i (t,x)=J24(x)Mt), 

k>0 

where {a k , k > 1} is an orthonormal basis in L2(M rf ; With v as in (|12.2U|) . the passive 
scalar equation ()12.4)) becomes 

(12.21) 9{t, x) = vA9(t, x) + f{t, x) - V0(t, x) ■ W{t, x), 

where W = W(t, x) is a d-dimensional space-time white noise and the ltd stochastic differ- 
ential is used. Previously, such equations have been studied using white noise approach in 
the space of Hida distributions |4l I4U| . A summary of the related results can be found in 
[H Section 4.3]. 

The Q-weighted Wiener chaos spaces allow us to state a result that is fully analogous to 
Theorem 112. II The proof is derived from Theorem I9.lt see j2H] for details. 

Theorem 12.6. Suppose that v > is a real number, each \a k (x)\ is a bounded mea- 
surable function, and the input data are deterministic and satisfy uq G L2(K rf ), / G 
L 2 ((0,r);ff 2 - 1 (K d )). 

Fix e > and let Q = {q^, k > 1} be a sequence so that, for all x, y G W 1 , 

Mv? - Yl lW k (x)al(x)y iyj > s\y\ 2 . 

k>\ 

Then, for every T > 0, there exits a unique w(H2(M. d ), H^ 1 ^ 1 )) Wiener Chaos solution 9 
of equation 

(12.22) d9(t, x) = (vA9(t, x) + f(t, x))dt - a k (x) ■ V9(t, x)dw k {t), 
The solution is an J-t-adapted process and satisfies 

ll lli 2iQ (W;L 2 ((O,T) ; H 2 1 (M<i))) + \\ d \\L 2iQ (W;C((0,T);L 2 (Rd))) 

< C(u,q,T) (|Nli 2(Rd) + \\f\\l 2ao ,Ty,H^(^))) ■ 

If maxj sup x \cr k (x)\ < C k , k > 1, then a possible choice of Q is 

q k = {5v) l/2 /(d2 k C k ), Q<5<2. 

If <rji.(x)<x[(:r) < Cfj < +oo, i,j = l,...,d,xE M d , then a possible choice of Q is 

q k = e{2u/{C (T d)) 112 , 0<e<l. 

13. Stochastic Navier-Stokes Equation 

In this section, we review the main facts about the stochastic Navier-Stokes equation and 
indicate how the Wiener Chaos approach can be used in the study of non-linear equations. 
Most of the results of this section come from the two papers [SHI and |31j . 

A priori, it is not clear in what sense the motion described by Kraichnan's velocity (see 
Section I12|) might fit into the paradigm of Newtonian mechanics. Accordingly, relating 
the Kraichnan velocity field v to classic fluid mechanics naturally leads to the question 
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whether we can compensate v (t, x) by a field u (t, x) that is more regular with respect to 
the time variable, so that there is a balance of momentum for the resulting field U (t, x) = 
u (t, x) + v (t, x) or, equivalently, that the motion of a fluid particle in the velocity field 
U (t, x) satisfies the Second Law of Newton. 

A positive answer to this question is given in jHS]) where it is shown that the equation for 
the smooth component u = (u 1 , . . . , u d ) of the velocity is given by 

du* = [vAu* - uWjU 1 - DiP + fi]dt 
(13.1) { +{g{- DiP k - D^u^) dw k , i = l,...,d, < t < T; 

divu = 0, u(0, x) = uo(x). 

where Wk, k > 1 are independent standard Wiener processes on a stochastic basis F, 
the functions a J k are given by (|12.3j) . the known functions f = (f 1 , . . . , f d ), gk = (<7^,), 
i = i, . . . ,d, k > 1 are, respectively, the drift and the diffusion components of the free force, 
and the unknown functions P, Pk are the drift and diffusion components of the pressure. 

Remark 13.1. It is useful to study equation M3.1\) for more general coefficients a\. So, in 
the future, a° k are not necessarily the same as in Section U^ 

We make the following assumptions: 
NS1 The functions at = ai(t,x) are deterministic and measurable, 



k>l \i=l ) 



< K, 



and there exists e > so that, for all y £ 



t G [0,T], x G R d . 
NS2 The functions /*, gi are non-random and 




i II 2 || i || 2 I 

i2((0,T);HJ 1 (M<*)) + 2^ ll^llt 2 ((0,T);i 2 (K d )) ) < °°- 



k>l 



Remark 13.2. In NS1, the derivatives Dio\ are understood as Schwartz distributions, but 
it is assumed that diva '■= ^2i = i diO~ l is a bounded I2 — valued function. Obviously, the latter 
assumption holds in the important case when Yli=i ^i a% = 0- 

Our next step is to use the divergence-free property of u to eliminate the pressure P and 
P from equation (|13.1|) . For that, we need the decomposition of ^(IR^M^) into potential 
and solenoidal components. 

Write 6(L 2 (R d ;lR d )) = {V G L 2 (M d ;lR d ) : div V = 0}. It is known (see e.g. H5J) that 

L 2 {R d ;R d ) = ®{L 2 (R d ;R d )) 6(L 2 (R d ;R d )), 
where 0(L2(M d ; R d )) is a Hilbert subspace orthogonal to &(L 2 
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The functions ©(V) and S(V) can be denned for V from any Sobolev space H^Mf 1 ; M. d ) and 
are usually referred to as the potential and the divergence free (or solenoidal), projections, 
respectively, of the vector field V. 

Now let u be a solution of equation (jl3.1|) . Since div u = 0, we have 

D^uAu* - uWjtf - D t P + f) = 0; D^D^u* +g\- A-Pfc) = 0, k > 1. 
As a result, 

DiP = 0(^An* - vPDjut + /*); D l P k = e^DjU* + g{), i = l,...,d, k > 1. 

So, instead of equation (|13,1[) . we can and will consider its equivalent form for the unknown 
vector u = (u 1 , . . . , u d ): 

(13.2) du = 6(vAu - vPDju + f)dt + 6(a J k DjU + g k )dw k , < t < T, 

with initial condition u|t=o = uo- 

Definition 13.3. An Tt-adapted random process u from the space Z^^ix [0, T]; H\ (R d ; M. d )) 
is called a solution of equation \13.$fy if 

(1) With probability one, the process u is weakly continuous in L2(M. d ;M. d ). 

(2) For every (p 6 Cq°(IR c( , M. d ), with div ip = there exists a measurable set SI' C SI so 
that, for all t £ [0,T], the equality 



(uV)(£) = K,<^)+ / {{vD jU \Drf){ S ) + (f\^){s))ds 

(1.3.3) 

[o-lDjU 1 + g\ip i )dw k {s) 

holds on SI'. In \V3. 3]) . (•,•) is the inner product in L2(J& d ) and (-,-,) is the duality 
between H^(R d ) and H^ l (W i ). 

The following existence and uniqueness result is proved in |31| . 

Theorem 13.4. In addition to NS1 and NS2, assume that the initial condition uo is non- 
random and belongs to L2(M d ; Then there exist a stochastic basis F = (SI, {J r t}t>o^) 
with the usual assumptions, a collection {w k ,k > 1} of independent standard Wiener pro- 
cesses on F, and a process u so that u is a solution of 113. and 



E [ sup ||u(s)||2 2(Rd;Rd) + / ||Vu(s)||| 2(Rd:Rd ^ ds ) < oo. 



If, in addition, d = 2, then the solution of 113.2(1 exists on any prescribed stochastic basis, is 
strongly continuous in t, is -adapted, and is unique, both path-wise and in distribution. 

When d > 3, existence of a strong solution as well as uniqueness (strong or weak) for 
equation (|13.2j) are important open problems. 

By the Cameron-Martin theorem, 

u(t,x) = u a (t,x)£ a . 

If the solution of (|13.2|) is J-^-adapted, then, using the ltd formula together with relation 
(|5.5[) for the time evolution of E(£ a | J^) and relation (j!2.17f) for the product of two elements 



WIENER CHAOS FOR STOCHASTIC EQUATIONS 57 

of the Cameron-Martin basis, we can derive the propagator system for coefficients u a [2D 
Theorem 3.2]: 



Theorem 13.5. In addition to NS1 and NS2, assume that uo £ £2(1^; R d ) and equation 
2]) has an -adapted solution u so that 



(13.4) supE||u||^ 2(Rd . Rd) (i) < 00. 
Then 

(13.5) u (t, x) = y^u a (t, x) £ a , 

and the Hermite- Fourier coefficients u a (t, x) are L2(M. d ; W d )-valued weakly continuous func- 
tions so that 

(13.6) sup \\ u a\\l 2{R d. Rd) (t) + Yl H Vu «llL(K d ;R rf x d )^) d * < °°- 
The functions u a (t, x) , a £ J , satisfy the (nonlinear) propagator 

^u Q = 6(Au a - Y *(Q!,/3i7) ( u 7' Vu (?) + J {|a|=o}f 

7,/3GA(a) 

(13.7) 

+ J2 \M ( V* V ) Ua ~^ + 7 {M=i}8*) ^ (t) ) , < t < T; 

U a | t= o = U /{| Q | = o}; 

recall that the numbers ^(ck, (3, 7) are defined in 112.1!$) . 

One of the questions in the theory of the Navier-Stokes equation is computation of the mean 
value u = Eu of the solution. The traditional approach relies on the Reynolds equation for 
the mean 



(13.8) d t u - i/Au + ( u, V) u = 0, 

which is not really an equation with respect to u. Decoupling (j!3.8|) has been an area of 
active research: Reynolds approximations, coupled equations for the moments, Gaussian 
closures, and so on (see e.g. [3^], |45j and the references therein) 

Another way to compute u (t, x) is to find the distribution of v (t, x) using the infinite- 
dimensional Kolmogorov equation associated with Q13.2|) . The complexity of this Kol- 
mogorov equation is prohibitive for any realistic application, at least for now. 

The propagator provides a third way: expressing the mean and other statistical moments 
of u in terms of u a . Indeed, by Cameron-Martin Theorem, 

Eu(i, x) = uo(t, x), 

Eu l (t,x)u (s,y) = ^<(i,x)<(s,y) 

If exist, the third- and fourth-order moments can be computed using ()12.15|) and (|12.16|) . 

The next theorem, proved in |31| , shows that the existence of a solution of the propagator 
(|13.7j) is not only necessary but, to some extent, sufficient for the global existence of a 
probabilistically strong solution of the stochastic Navier-Stokes equation ()13.2|) . 
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Theorem 13.6. Let NS1 and NS2 hold and u G L 2 (R d ;R d )- Assume that the propagator 
\13. ?D has a solution {u a (t, x) , a € J} on the interval (0, T] so that, for every a, the 
process u a is weakly continuous in L2(M. d ;M. d ) and the inequality 

(13.9) sup 1 1 "a 1 1 1 2 (]£<*;»<*)(*) + / Yl H Vu «lli 2 (IR d ;IR d x d )^) dt < 00 
holds. If the process 

(13.10) U(t,x) := u Q (i,x)£ a 

is -adapted, then it is a solution of \V3. 
27ie process tj satisfies 



d) ds J < oo 



E^up||U( S )||2 2(Rd;Rd) + ^ ||VU( S )||i 2(Rd 
and, /or every v E L 2 (M d ;M d ), E (tj, v) is a continuous function oft. 



Since U is constructed on a prescribed stochastic basis and over a prescribed time interval 
[0, T] , this solution of (|13.2|) is strong in the probabilistic sense and is global in time. 
Being true in any space dimension d, Theorem 113.61 suggests another possible way to study 
equation (J13.2|) when d > 3. Unlike the propagator for the linear equation, the system 
()13.7|) is not lower-triangular and not solvable by induction, so that analysis of ()13.7j) is an 
open problem. 



14. First-Order Ito Equations 

The objective of this section is to study equation 

(14.1) du(t,x) = u x (t,x)dw(t), t > 0, x G R, 
and its analog for x 6 M. d . 

Equation (|14.1I) was first encountered in Example l6.8l see also j^j . With a non-random initial 
condition u(0,x) = <p(x), direct computations show that, if exists, the Fourier transform 
u = u(t, y) of the solution must satisfy 

(14.2) du(t,y) = V^lyu(t,y)dw(t), or u(t,y) = ^e^^+^K 

The last equality shows that the properties of the solution essentially depend on the initial 
condition, and, in general, the solution is not in Z/2(W). 

The S-transformed equation, v t = h(t)v x , has a unique solution 

v (t, x) = ip (^x + J h(s)ds^j , h(t) = hjiriiit). 

The results of Section |3] imply that a white noise solution of the equation can exist only if 
(p is a real analytic function. On the other hand, if 92 is infinitely differentiable, then, by 
Theorem 18.41 the Wiener Chaos solution exists and can be recovered from v. 
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Theorem 14.1. Assume that the initial condition tp belongs to the Schwarz space S = <S(R) 
of tempered distributions. Then there exists a generalized random process u = u(t, x), t > 0, 
x G R, so that, for every 7 G R and T > 0, i/ie process u is the unique u;(i?2 (R), fl^" 
Wiener Chaos solution of equation HXTp. 



Proof. The propagator for Q14.1JI is 

(14.3) u a (t,x) = ip(x)I{\a\ = 0) + / V" v ^(it a - ( j)(s,x)) :c m i (s)ds. 

Even though Theorem 16.41 is not applicable, the system can be solved by induction if <p 
is sufficiently smooth. Denote by C ip (k), k > 0, the square of the -L 2 (R) norm of the A: th 
derivative of (p: 



-00 



(14.4) C v (fc) = / |<^(x)| 2 d:r. 

J —00 

By Corollary 16.61 for every > and n > 0, 

(14.5) Ell(^ ) )*llW*) = < * C ^ + fc) - 

|a|=fc 

The statement of the theorem now follows. □ 

Remark 14.2. Once interpreted in a suitable sense, the Wiener Chaos solution of \14-l\) 
is -adapted and does not depend on the choice of the Cameron- Martin basis in ^(W). 
Indeed, choose the wight sequence so that 

1 



r 



,2 



" " l + C v (\a\ 



By J773| ) ; we have u G KL 2 (W; L 2 (R)). 
ATexf, define 



1 sinfiVa;) 

^at(x) - 



7T X 

Direct computations show that the Fourier transform of tp^ is supported in [— N, N] and 
f-^^N(x)dx = 1. Consider equation \14-1\) with initial condition 

Pn(x) = / p(x -y)^ N {y)dy. 

By [14-ty , this equation has a unique solution so that u^(t, •) G L 2 (W; iT 2 (R)), t > 0, 
7 G R. Relation \14-£\) and f/ie definition of un imply 

N^oo ^ " Ma ~ UjV ' a llL(R)(*) = °' * - °' fc - °' 
|«|=fc 

so t/iat, fry t/ie Lebesgue dominated convergence theorem, 

lirn^ \\u - Mjv|||L 2 (W;L 2 (R))0) = °> * > 0. 

In other words, the solution of the propagator j!4-'A ) corresponding to any basis m in 
Li2((0,T)) is a limit in IZL2 (W; L 2 (R)) of the sequence {un, N > 1} of J 7 ^ -adapted pro- 
cesses. 



The properties of the Wiener Chaos solution of (|14.1|) depend on the growth rate of the 
numbers C v {n). In particular, 
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• If C v {n) < C n (niy, C>0, < 7 < 1, then 

u<EL 2 (W; L 2 ((0, T); flJ(R))) for all T > and every n > 0. 

• If C<p(n) < C"n!, C>0, then 

- for every n > 0, there is a T > so that u G L 2 (W; L a ((0, T); J3J(R))). In 
other words, the square-integrable solution exists only for sufficiently small T. 

— for every n > and every T > 0, there exists a number 5 G (0, 1) so that 
u G L 2iQ (W; L 2 ((0, T); H r 2 l (R))) with Q = («5, 5,5,...). 

• If the numbers CL(«) grow as C n (n!) 1+P , p > 0, then, for every T > 0, there exists 
a number 7 > so that 

u G (S)-p- 7 (L 2 (W);L 2 ((0,r);fly(R))). If p > 0, then this solution does not 
belong to any L 2)Q (W; L 2 ((0, T); i?£(R))). If p > 1, then this solution does not 
have an S-transform. 

• If the numbers C^(n) grow faster than C n {n\) b for any b, C > 0, then the Wiener 
Chaos solution of (|14.1(l does not belong to any 

(S)_ p ,_ 7 (L 2 ((0,T);# 2 "(R))), p, 7 > 0, or L 2 , Q (W;L 2 ((0,T);flJ(R))). 
To construct a function <p with the required rate of growth of C^n), consider 

/>oo 

ip(x) = / cos(xy)e~ 9 ^dy, 
Jo 

where g is a suitable positive, unbounded, even function. Note that, up to a multi- 
plicative constant, the Fourier transform of ip is e~ 9 ^ y \ and so C^n) grows with n as 
f +C °\y\ 2n e~ 2 ^dy. 

A more general first-order equation can be considered: 

(14.6) du(t,x) = a ik (t,x)Diu(t,x)dw k (t), t > 0, x G M d . 

Theorem 14.3. Assume that in equation \14-b\) the initial condition n(0, x) belongs to 
S(R d ) and each is infinitely differentiable with respect to x so that &\xpr tx \ \D n aik(t, x)\ < 
Cik{n), n > 0. Then there exists a generalized random process u = u(t, x), t > 0, x G M. d , so 
that, for every 7 G R and T > 0, the process u is the unique w(H2(M. d ), (M. d )) Wiener 
Chaos solution of equation {l^.lty . 

Proof. The arguments are identical to the proof of Theorem 114.11 □ 

Note that the S-transformed equation (|14.6j) is v% = h^Gi^DiV and has a unique solution if 
each aik is a Lipschitz continuous function of x. Still, without additional smoothness, it is 
impossible to relate this solution to any generalized random process. 
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